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1.  EXECUTIVE  SUMMARY 


The  research  investigations  contained  in  this  report  are  concerned  with 
thermoelastic  behavior  of  cracked  composite  materials  and  with  analysis  of  damage 
formation  in  the  forms  of  cracks  and  debonding  of  the  interface  between  the  constituents. 
In  the  first  case  the  damage  is  geometrically  defined  ,  for  example  in  the  form  of  crack 
distributions  or  degraded  interface,  and  the  goal  is  to  evaluate  the  effect  of  such  damage 
on  the  stiflfiiess  and  the  thermal  expansion  coefficients  (TEC)  of  the  composite  material.  In 
the  second  case  the  goal  is  to  predict  the  crack  density  which  is  produced  within  a  certmn 
composite  by  specified  load  and  temperature  inputs  and  to  evaluate  load/temperature 
inputs  which  produce  complete  debonding  of  the  constituents. 

There  are  five  parts  as  listed  in  the  table  of  contents.  The  first  three  are  concerned 
with  evaluation  of  thermoelastic  properties  and  the  last  two  are  concerned  with  the 
evaluation  of  damage .  These  investigations  will  now  be  discussed  sequentially. 

Cracked  Laminates  with  Imperfect  Interiaminar  Interface 

This  work  deals  with  analysis  of  laminates  with  imperfect  interlaminar  interface. 
Such  an  imperfect  interface  models  the  case  when  there  is  a  thin  layer  in-between  the  plies 
of  the  laminate,  a  situation  encountered  when  there  are  oxidation  protection  layers,  or 
when  the  interface  has  suffered  damage  for  whatever  reasons.  The  mathematical 
description  of  such  imperfection  is  a  discontinuity  of  displacement  which  is  related  to  the 
interface  traction.  In  our  case  the  relation  between  interface  displacement  discontinuity 
and  traction  is  linear.  Such  a  relation  is  caUed  an  imperfect  interface  condition. 

It  should  be  noted  that  in  the  common  case  of  a  symmetric  laminate  subjected  to 
membrane  loads  in  it’s  plane  the  imperfect  interface  has  no  effect  ,  neither  on  internal 
stresses  nor  on  the  laminate  thermoplastic  properties.  But  the  present  analysis  shows  that 
the  situation  is  entirely  different  for  a  laminate  which  has  intralaminar  cracks.  The  method 
of  analysis  is  an  extension  of  a  previous  one  which  dealt  with  analysis  of  cracked  laminates 
with  perfect  interface  [1].  This  analysis  was  based  on  an  application  of  the  classical 
variational  principle  of  minimum  complementary  energy.  This  minimum  principle  has  been 
extended  to  the  case  of  imperfect  interface  in  [2]  and  this  extended  principle  has  been 
applied  for  the  present  situation. 

The  analysis  shows  that  for  ceramic  composite  cross-ply  laminates  interface 
imperfection  has  only  a  small  effect  on  the  in-plane  Young’s  modulus  and  TEC  of  the 
composite  but  has  a  significant  effect  on  the  laminate  shear  modulus  and  the  interlaminar 
shear  and  normal  stresses. 

Differential  Scheme  for  Effective  Thermoelastic  Properties  of  Cracked  Composite 
Materials. 

The  most  common  damage  in  composite  materials  is  in  the  form  of  microcracks. 
The  research  described  in  paragraph  1  above  was  concerned  with  microcracks  through 
laminate  ply  thickness.  Another  very  important  problem  is  the  effect  of  matrix  microcracks 
in  a  unidirectional  fiber  composite  or  in  a  particulate  composite.  The  size  of  such  micro 


cracks  is  of  the  order  of  fiber  or  particle  diameters.  Essentially,  such  a  cracked  composite 
is  a  three  phase  material  consisting  of  matrix,  reinforcement  (fibers  or  particles)  and 
cracks,  the  goal  of  the  research  is  to  evaluate  the  effect  of  such  microcracks  on 
thermoelastic  properties,  that  is  stifi&iess  and  TEC.  There  exists  a  veiy  extensive  literatare 
on  evaluation  of  thermoelastic  properties  of  reinforced  materials  and  of  cracked  materials. 
The  difficulty  of  the  problem  is  demonstrated  by  the  fact  that  there  exists  almost  no 
literature  on  evaluation  of  properties  of  cracked  composites  of  the  kind  described. 

The  problem  can  be  solved  exactly  only  in  the  case  Avhen  the  composite  is 
described  by  a  matrix  containing  a  small  amount  of  non  interacting  fibers/particles  and 
cracks.  But  such  a  problem  is  only  of  academic  interest.  It  is  therefore  necessary  to  resort 
to  approximate  methods.  Two  well  known  methods  of  approximation  are  the  self 
consistent  scheme  (SCS)  and  the  differential  scheme  (DS).  The  SCS  has  actually  been 
applied  to  the  problem  under  consideration  ,however  ,  it  is  well  known  that  the  SCS 
results  for  cracked  homogeneous  materials  gives  unreliable  results  in  that  it  seriously 
underestimates  the  stiflfiiess.  It  can  therefore  not  be  expected  to  yield  reliable  results  for 
cracked  composites.  On  the  other  hand  the  DS  results  for  cracked  homogeneous  materials 
are  much  more  acceptable  and  also  agree  vdth  experiments  performed.  It  appeared 
therefore  reasonable  to  apply  this  method  for  the  present  problem. 

The  basic  premise  of  the  DS  is  that  when  a  small  number  of  inhomogeneities 
(fibers,particles,cracks)  are  added  within  a  composite  then  their  incremental  effect  is  given 
by  dilute  concentration  theory  vrith  respect  to  the  current  properties  of  the  composite.  The 
DS  becomes  a  quite  accurate  method  when  the  added  inhomogeneities  are  by  an  order  of 
magnitude  larger  than  the  ones  previously  added. 

The  present  investigation  has  resulted  in  evaluation  of  effective  elastic  moduli  and 
TEC  of  a  composite  containing  particles  or  fibers  and  matrix  cracks.  To  the  best  of  our 
knowledge  this  the  first  time  that  the  TEC  of  this  kind  of  composite  has  been  analyzed. 

Macro-Residual  Strains  due  to  Cyclic  Loading  of  Composites 

It  is  well  known  that  when  ceramic  composites  are  subjected  to  load  cycles,  at 
room  temperature,  macro-residual  strains  remain.  Such  residual  strains  are  commonly 
observed  in  metals  which  have  been  loaded  into  a  plastic  state  and  then  unloaded. 
However,  ceramic  composites  are  elastic-brittle  and  thus  the  phenomenon  requires 
explanation  .  It  is  shown  that  such  residual  strains  are  formed  because  of  the  residual 
stress  state  in  the  composite  which  is  due  to  the  fact  that  the  composite  is  stress  free  at  it’s 
high  formation  temperature  and  thus  develops  residual  stress  due  to  cooldown.  If  the 
loading  part  of  the  load  cycle  produces  internal  cracks  then  the  residual  stress  state 
changes.  The  stresses  due  to  the  load  disappear  on  unloading  but  the  composite  is  now  in 
a  state  of  different  residual  stress  and  strain  and  therefore  the  average  or  macro-strain  has 
changed  after  unloading.  The  difference  between  the  average  strain  before  loading  and 
after  unloading  is  the  measured  macro-residual  strain. 

The  fundamental  novel  result  found  in  this  investigation  is  that  the  macro-residual 
strain  can  be  expressed  in  terms  of  the  difference  between  the  effective  TEC  after  and 
before  the  loading  cycle.  Furthermore,  this  explains  the  physical  phenomenon  observed 
that  after  many  load  cycles  the  residual  strain  will  achieve  a  limit  value  which  cannot  be 


increased  by  additional  load  cycles..  The  reason  for  this  is  that  crack  density  cannot  grow 
indefinitely  and  at  some  time  crack  saturation  will  take  place,  which  implies  that  no 
additional  cracks  can  be  introduced  into  the  material.  Now  the  effective  TEC  are  functions 
of  the  crack  geometry  and  thus  at  crack  saturation  the  effective  TEC  assumes  a  limit 
value.  Therefore,  the  macro-residual  str^  also  assumes  a  limit  value  since  it  is  uniquely 
determined  by  the  effective  TEC. 

Finite  Thermoelastic  Fracture  Criterion  with  Application  to  Laminate  Cracking 
Analysis 

The  origin  of  classical  fracture  mechanics  is  wth  a  differential  energy  balance.  The 
question  posed  is  :  what  is  the  load  which  pennits  differential  extension  of  a  crack  without 
load  increase  ?  The  answer  is  given  by  an  energy  balance  in  which  the  change  in  external 
work  due  to  crack  extension  is  equated  to  the  internal  energy  change  plus  the  energy 
needed  to  open  up  the  crack.  The  concern  is  with  single  cracks  or  perhaps  with  a  number 
of  cracks  at  well  known  locations  wthin  a  homogeneous  material. 

In  composite  materials  in  general,  and  in  fiber  composite  laminates  in  particular, 
the  formation  of  damage  consists  of  the  accumulation  of  a  multitude  of  interacting  cracks. 
It  is  impossible  to  follow  the  differential  extensions  of  such  cracks  which  appear  suddenly 
and  at  unspecified  locations.  Instead  we  have  what  may  be  called  fracture  events  which 
may  be  assumed  to  occur  spontaneously.  Therefore  a  different  approach  has  been 
developed  which  defines  a  critical  load  as  one  which  permits  two  simultaneous  different 
crack  geometry’s  which  differ  by  finite  fi-acture  surface.  Therefore  the  approach  may  be 
termed  finite  fracture  mechanics. 

Another  consideration  which  is  peculiar  to  composites  is  that  there  are  residual 
stresses  which  are  a  result  of  cooldown  after  manufacturing  at  high  temperature  and 
constituent  TEC  mismatch.  Thus  load  produces  stresses  which  must  be  added  to  the 
residual  stresses. 

We  have  constructed  a  criterion  for  evaluation  of  the  critical  load  which  permits 
the  spontaneous  appearance  of  finite  new  fracture  surface.  This  criterion  is  expressed 
solely  in  terms  of  the  stress  fields  before  and  after  the  appearance  of  the  new  fracture 
surface.  Moreover,  it  has  been  shown  that  the  critical  energy  release  which  is  needed  to 
produce  the  new  fracture  surface  can  be  bounded  from  above  by  use  of  admissible  stress 
fields  instead  of  the  actual  stress  fields. 

This  theory  has  been  applied  to  prediction  of  intralaminar  crack  densities  in 
laminates  which  are  cooled  down  from  a  stress-free  temperature  and  are  subsequently 
subjected  to  mechanical  load.  For  this  purpose  a  novel  variational  method  for  analysis  of 
thermoplastic  stresses  in  laminates  has  been  constructed.  The  loads  which  initiate 
intralaminar  cracking  have  been  evaluate  for  various  cross-ply  laminates  .Good  agreement 
with  experimental  results  has  been  obtained. 

Interface  Debond  Analysis 

A  possible  failure  mechnism  in  composite  materials  is  debonding  of  the  interface 
between  the  constituents.  If  this  happens  the  reinforcing  fibers  or  inclusions  are  essentially 


replaced  by  voids  and  the  composite  material  is  replaced  by  a  porous  material.  A  porous 
material  is  weaker  than  the  matrix  and  cannot  serve  as  a  structural  material. 

The  debonding  process  consists  of  the  growth  of  interface  cracks  which  separate 
the  constituents.  Analysis  of  such  a  process  is  an  impossible  task  for  even  the  partial 
debonding  of  a  single  inclusion  isolated  in  an  infinite  matrix  is  a  very  formidable  problem. 
Here  we  consider  a  different  problem  :  what  is  the  load/  temperature  combination  which 
produces  debonding  of  the  entire  interface  ?  This  problem  can  be  treated  in  relatively 
simple  fashion  by  considering  it  in  terms  of  the  finite  fi-acture  criterion  theory  discussed  in 
par.  4,  above.  It  is  assumed  that  the  debonding  of  the  interface  is  a  spontaneous  fi-acture 
event.  Then  the  the  two  comparison  stress  states  are  the  stress  in  the  composite  for  perfect 
bond  and  the  stress  for  complete  debonding,  i.e.  a  porous  material.  These  stress  fields  can 
be  evaluated  for  suitable  models  of  a  composite.  In  the  present  study  we  have  chosen  the 
model  known  by  the  name  Generalized  Self  Consistent  Scheme. 

The  an^ysis  produces  a  temperature  dependent  debond  failure  surface  which  is 
expressed  in  terms  of  applied  average  stress  and  temperature.  The  validity  of  the  failure 
surface  is  limited  by  the  need  to  assure  that  after  debonding  particles  do  not  penetrate  the 
matrix.  The  theory  has  been  applied  to  evaluation  of  debond  failure  surfaces  of  a 
unidirectional  fiber  composite  which  is  loaded  transversely  with  temperature  change. 
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CRACKED  LAMINATES  WITH  IMPERFECT 
INTERLAMINAR  INTERFACE 


Z.  Hashin 


Cracked  Laminates  with  Imperfect  Interlaminar 
Interface 


1.  Introduction 

The  present  work  is  concerned  with  the  effect  of  intralaminar  crack  (IC) 
accumulation  on  the  thermomechanical  properties  of  fiber  composite  lam¬ 
inates  and  the  resulting  internal  stress  distributions.  Such  cracks  develop 
in  the  matrix  along  fibers  due  to  load  or  temperature  change.  They  are 
thus  parallel  crack  distributions  within  the  layers  which  propagate  very 
rapidly  until  the  laminate  edges.  Therefore, the  formation  of  a  typical  IC  is 
not  viewed  as  a  crack  propagation  phenomenon  but  as  a  fracture  event 
which  occurs  instantaneously.  Thus,  the  concern  is  with  a  laminate  which 
contains  IC  distributions  which  are  quantitatively  described  by  crack  den¬ 
sity,  the  number  of  IC  per  unit  length.  The  problems  are  then  to  determine 
deterioration  of  thermoelastic  peoperties  in  terms  of  crack  density,  lami¬ 
nate  internal  geometry  and  ply  properties,  internal  stresses  resulting  from 
crack  accumulation  and  their  relation  to  failure  mechanisms,  and  more  am¬ 
bitiously  -  to  predict  crack  density  due  to  load  or  temperature. 
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The  problems  outlined  have  been  the  subject  of  a  large  number  of  re¬ 
search  papers  over  the  last  15  years.  There  are  two  major  approaches  : 
the  first  may  be  termed  the  micromechanics  approach  and  the  second,  the 
continuum  damage  approach.  In  the  first  approach  it  is  attempted  to  car¬ 
ry  out  analysis  recognizing  the  cracks  as  defects  on  which  the  tractions 
must  vanish.  The  advantage  of  this  approach  is  physical  realism  and  in¬ 
formation  about  internal  (micro)  stresses,  which  is  important  for  fmlure 
considerations.  The  disadvantage  is  analytical  difficulty  and  for  this  reason 
the  micromechanics  approach  has  to  date  been  confined  to  cross-plies. 

In  the  second  approach  effect  of  IC  on  a  layer  is  modeled  by  an  abstract 
damage  function  whose  form  is  not  unique  and  which  invariably  contains 
unknown  coefficients.  The  disadvantage  is  that  such  coefficients  must  be 
backed  out  from  experiment  on  the  laminate  and  it  is  not  dear  whether 
such  coefficients  qualify  as  ply  material  parameters  or  are  fitting  parameters 
which  change  from  laminate  to  laminate.  The  advantage  of  the  approach 
is  that  it  can  be  applied  to  practical  laminates,  more  complicated  than 
cross-plies. 

The  present  work  is  concerned  with  the  micromechanics  approach  for 
cross-ply  laminates.  Review  of  the  voluminous  literature  is  not  within  the 
present  scope.  It  is  recalled  that  initial  analytical  efforts  were  based  on 
the  shear-lag  approximation  e.g.  Reifsnider  and  Jamison  (1982),Laws  and 
Dvorak  (1988).  This  method  requires  the  determination  of  a  so-called  shear 
lag  parameter  on  the  basis  of  the  fracture  toughness  of  the  ply  materials. 
Analysis  in  terms  of  a  displacement  formulation  represented,  arbitrarily, 
by  hyperbolic  functions  was  given  by  Tsai  et  al.  (1990).  Work  of  similar 
nature  with  the  choice  of  different  form  displacement  functions  has  been 
done  by  Lee  et  al.  (1990).  A  variational  method  based  on  the  prindple  of 
TtiiniTmiTn  complementary  energy  has  been  developed  by  Hashin  (1985)  with 
application  to  stiffness  reduction  and  stress  analysis  of  cross-ply  laminates 
with  one  layer  cracked.  This  has  been  extended  to  the  case  of  all  layers 
cracked  in  Hashin  (1987)  and  to  evaluation  of  thermal  expansion  coefficients 
in  Hashin  (1988).  The  only  assumption  made  in  the  variational  analysis  is 
that-in  plane  stresses  in  the  ply  are  constant  over  the  thickness.  Analysis 
based  on  similar  assumptions  has  been  given  by  McCartney  (1992).  Analysis 
for  more  general  in-plane  stresses  has  been  given  by  Varna  and  Berglund 
(1994).  Naim  et  al.  (1993)  have  successfully  used  the  variational  analysis 
for  prediction  of  crack  density  resulting  from  in  plane  loading  of  cross-ply 
laminates,  , 

The  purpose  of  the  work  presented  here  is  to  extend  the  variational 
analysis  to  the  evaluation  of  thermoelastic  properties  and  internal  stresses 
of  cross-ply  laminates  when  there  is  imperfect  interlaminar  bond  between 
the  layers. 
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2.  Thermo  elastic  Extremum  Principle  for  Imperfect  Interface 

Perfect  interface  between  two  solid  constituents  implies  continuity  of  trac¬ 
tion  and  displacement  vectors  at  the  interface.  When  the  interface  displace¬ 
ment  vector  is  discontinous,  while  the  traction  vector  remains  continuous 
for  reasons  of  equilibrium,  the  interface  is  called  imperfect.  Let  the  dis¬ 
placement  jump  at  interface  S12  be  denoted 


[u]  =  —  u^. 

(1) 

Then  the  simplest  imperfect  interface  condition  is 

Tn  =  DnW 

r,  =  Ds  Kl 

(2) 

Tt  =  Dt  N, 

where  n,s,<  are  normal  and  tangential  components  of  the  interface  normal 
n,  assumed  here  as  pointing  into  phase  2,  and  Dnt  spring  constant 

type  interface  parameters.  With  respect  to  a  fixed  cartesian  coordinate 
system,  (2)  assumes  the  forms 

T  =  D.lu]  [u]  =  R.T  R==D-S  (3) 

where  the  Cartesian  components  of  D  and  its  inverse  R  now  vary  along  the 
interface.  It  has  been  shown  in  Hashin  (1990)  that  the  effect  of  a  thin  and 
very  compliant  interphase  between  constituents  can  actually  be  expressed 
in  the  form  (2)  and  that  the  interface  parameters  can  be  expressed  in  terms 
of  interphase  thickness  and  stiffiaess. 

In  the  variational  analysis  to  be  employed  here  the  generalization  of 
the  extremum  principle  of  minimum  complementary  energy  for  imperfect 
interface  conditions  will  be  needed,  Hashin  (1992).  This  will  here  be  further 
generalized  to  the  thermoelastic  case  and  will  be  stated  for  the  case  when 
tractions  are  prescribed  over  the  entire  external  surface  S.  Let  <r  be  the 
actual  stress  field  and  a  an  admissible  stress  field  for  a  body  with  surface 
load  T  (5)  and  imperfect  interface  512.  Define 


W  = 

.W  =  I? 


S 

S 


(4) 


where  S  is  the  compliance  tensor.  Here  W  is  the  stress  energy  density  while 
W  has  no  physical  meaning.  Next  define  the  functionals 


U  =  /j,[W-l-«.^^-Cp(5V20o)]dF-l-|4^T:R:Td5 
U  =  /y[w  +  a.?5-Cp(0V2flo)]dt^+l4,T:R:Td5 
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Here  «  is  the  thermal  expansion  tensor,  Cp  the  specific  heat  at  constant 
pressure  and  0  is  the  (known)  temperature  relative  to  a  reference  tem¬ 
perature  0o.  Then  the  thermoelastic  principle  of  minimum  complementary 
energy  is  expressed  by  the  inequality 

u>ir,  (®) 

equality  occurring  if,  and  only  if,  a  =  a  . 

For  composite  materials  applications  it  is  of  importance  to  consider  the 
case  of  constant  temperature  and  so-called  homogeneous  traction  boundary 
conditions  which  are  defined  as 

T(S)  ==  <r''.n(5)  (?) 

Where  o*®  is  a  constant  stress  tensor.  Then  cr®  is  the  average  stress  tensor 
and  it  can  be  shown  that  the  first  of  (5)  is,  rigorously 

H  i  [o-"  :  S*  :  O'"  -H  a*,  -  c*  V  (8) 

2  i 

where  S*  ,a*  and  Cp  are  the  eflfective  elastic  compliance  tensor,  thermal 
expansion  tensor  and  specific  heat,  respectively. 

In  the  following  the  variational  principle  will  be  exploited  to  analyze 
approximately  thermo-elastic  properties  and  internal  stresses  in  cracked 

laminates. 

S.  Cross-Ply  Laminates  with  One  Ply  Family  Cracked 

The  case  to  be  considered  here  is  a  [0®  ,  90®].  laminate  in  whi  A  either  the  0® 
or  the  90®plies  are  cracked,  fig.  la.  The  variational  method  will  be  employed 
to  obtain  strict  lower  bounds  for  the  effective  Yofing’s  modulus  F*  and  the 
effective  shear  modulus  G*^  and  approximations  for  the  effective  thermal 
expansion  coefficient  a*  and  internal  stresses,  for  the  case  of  imperfect 
interlaminar  interface  as  defined  by  a  damaged  interphase  between  plies. 
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Figure  1.  Cracked  laminate 


3.1.  EFFECTIVE  YOUNG’S  MODULUS 

Let  it  be  assumed  that  the  laminate  is  subjected  to  a  constant  tensile 
membrane  force  and  let  the  stresses  in  the  layers  of  the  uncracked 
laminate  be  denoted  Oi  and  02,  respectivdy,  where  from  now  on  the  label 
1  refers  to  the  90“  ply  and  2  refers  to  the  0“  ply.  As  is  weU  known,  these 
stresses  are  constant  throughout  the  layers  .  The  actual  stress  state  in  the 
cracked  laminate  is  described  by  generalized  plane  strain  in  reference  to  the 
y  axis,  and  therefore  all  stresses  are  functions  of  x,z  only.  It  is  at  present 
assumed  that  only  the  90“  ply  is  cracked.  The  admissible  stress  state  in  the 
cracked  laminate  udll  be  Constructed  on  the  basis  of  the  simplification  that 
the  Oxx  stresses  are  functions  of  x  and  not  of  z.  Thus  these  stresses  may  be 
written  in  the  form 

=  cTiIl  -  ^i(x)] 

(9) 


=  t^2[l  —  ^2(3;)] 


where  and  ^2  are  unknown  functions.  These  functions  are,  however, 
related  since  the  stress  pairs  ^2  and  (9)  are  each  in  equilibrium  with 
the  same  Therefore 


+  <^2^24^2  —  0 

Consider  a  typical  region  between  two  adjacent  cracks  at  distance  2a,  fig.lb. 
It  is  emphasized  that  the  cracks  do  not  have  to  be  equidistant.  For  example, 
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the  interdistance  2a  could  be  a  random  variable.  For  a  typical  intercra.ck 
redon  as  shown  in  fig.lb  the  admissible  stress  field  is  constructed  by  in¬ 
tegration  of  the  two  dimensional  equations  of  equilibrium  in  the  xz  plane 
with  the  stresses  (9).  The  integration  produces  residual  unknown  function- 
E  which  are  determined  by  satisfaction  of  traction  continuity  conditions 
at  layer  interfaces  and  zero  traction  conditions  on  the  external  lanunate 
surface.  The  remaining  admissible  stresses  are  then 


aif  =  ia^/\)4>'ix)(h-z) 

aiv  = 

where  (f>  =  4>i  ,  \  =  hlh  and  prime  superscript,  denotes  x  differentiation. 
These  stresses  have  already  been  given  in  Hashin  (1985).  On  the  crack 
surfaces  and  aiV  must  vanish  and  therefore 

^(±a)  =  1  ~  ® 

The  aggregate  of  the  stress  fields  (9-11)  for  all  intercrack  regions  are  the 
admissible  stress  field.  The  stress  energy  densities  for  the  layers  1  and  2 
are: 

2Wi  =  fEr  -  2a^}a^z)vTl Er  +  o^z)  /J^r  +  tr^z 

2W2  =  a^lEx  —  2a^}a^)vAlEA-\-o^z  jEr-^-ciz  JGa 

where  the  elastic  ply  properties  in  (12)  are :  Ea,  va  -  axial  Young’s  modulus 
and  associated  Poisson’s  ratio  ;  Er,  VT  -  transverse  Young’s  modulus  and 
associated  Poisson’s  ratio  ;  Gai  Gt  -  axial  and  transverse  shear  moduli. 

For  reasons  of  symmetry  it  is  sufficient  to  evaluate  the  complementary 
energy  functional  for  the  regions  — x  <  Umi  0<z</i. 

Then  for  such  a  region  and  for  an  isothermal  state,  the  second  of  (5)  assumes 

the  form 

fe™  =  S-Z.  C  widz<b + cr.  /;  mdzdx-h 

I  f-Z  ^1)/ 

where  I?„  and  Ds  are  normal  and  shear  interface  parameters,  and  for  the 
entire  laminate 

(14) 
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The  actual  complementary  energy  is  given  for  the  present  case  by 

Uc  =  o®  =  N^j2h  (15) 


where  V  is  the  volume  of  the  entire  cracked  laminate  of  thickness  1  in  y 
direction.  Introduction  of  the  admissible  stress  field  with  functions  4pm  lor 
each  region  into  (14)  results,  according  to  the  principle  of  minimum  corn- 
plementary  energy,  in  an  upper  bound  on  (15).  The  lowest  upper  bound  is 
obtained  by  minimizing  the  resulting  functional  with  respect  to  the  func¬ 
tions  <pm’  This  is  a  standard  problem  in  the  calculus  of  variations  resulting 
in  Euler  equations  and  boundary  conditions  for  the  minimizing  functions 
wMch  have  the  form 


+  P 


de 


+ 


=  0, 


^m(±/»m)  =  1> 


d<pm 


=  0 


(16) 


where 


f  =  x/ti  Pm  =  O-m/h  P  =  (<^02  "  Cn)IC22  9  -  C'oo/C'22 

Coo  =  iIJSt  +  iI>^Ea  Oo2  =  :^(a+|)-^a 

C22  =  (A  +  l)(3A2-H2A  +  8)/60Er+AV4I?„<i  (17) 

Cn  =  |(l/Gr  +  l/AG'A)+l/i?s<i 

Evaluation  of  the  complementary  energy  functional  in  terms  of  the  func¬ 
tions  <pm  as  was  done  in  Hashin(1985),  introduction  of  the  result  into  the 
complementary  inequality  (6)  with  (15)  as  the  actual  complementary  ener- 
gy  gives  the  result 


1 /;?*  <  =  (18) 
l/^*-l/^-  +  (J  A+1  (p)  ’  ^  ^ 


where  the  brackets  denote  average  with  respect  to  the  random  variable  Umj 
half  of  the  intercrack  spacing.  The  form  of  x  depends  on  the  nature  of  the 
roots  of  the  characteristic  equation  of  (16).  When  these  roots  are  of  the 
form  ±(a  -h  ifi),  where  i  =  ,  then  the  solution  of  (16)  is  of  the  form 

<pm  =  AmCosh(a()  cos(/3^)  -1-  BmSinh{a^)  sin(^^)  (19) 

where  the  constants  are  determined  by  the  boundciry  conditions  in  (16). 
Then  the  associated  Xm  is 


Xm  =  2a/3(a^  +  ^^) 


Cosh(2apm)  -  cos(2Ppm) 
a  sin(2/3pm)  +  fiSinh(2apm) 


(20) 
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If  the  roots  are  real,  thus  of  the  form  ±a,  ±/3  ,theii 
—  A.fnC  osh(^oc^'^  +  Bm.Goshifi^') 


(21) 


3.2.  IMPERFECT  INTERFACE  MODEL 

A  physical  interpretation  of  the  interface  parameters  I>„  and  D,  has  been 
given  in  Hashin  (1990).  If  there  is  a  thin  elastic  isotropic  interphase  of 
thickness  U  between  the  phases,  then 

D^  =  {Ki  +  \Gi)lti  D,  =  Dt=^Gilti  (22) 

where  Ki  and  Gi  are  the  bulk  and  shear  moduli  of  the  interphase.  It  is 
easUy  shown  that  if  the  interphase  is  orthotropic,  with  material  axes  n,s 
and  t  ,  (22)  becomes 

x>„  =  Cnn/ti  D,  =  Gn,/t.  Dt  =  G„t/ti  (23) 

where  Cnn  is  the  normal  stiffiiess  and  Gns  are  shear  moduli.  The 

first  of  (22,23)  is  strictly  valid  only  when  the  interphase  elastic  moduli  mc 
much  smaller  than  those  of  the  constituents,  but  there  is  no  such  restric¬ 
tion  with  respect  to  D,  and  Df  If  the  interphase  moduli  are  of  the  ord^ 
of  constituent  moduli  then  the  thin  interphase  effect  is  negligible  and  is 
equivalent  to  a  perfect  interface  with  displacement  continuity.  Consider  a 
thin  interphase  in-between  the  layers  of  a  cross-ply,  fig.2.  A  relevant  ex¬ 
ample  is  an  oxidation  protection  layer  between  the  laminae  of  a  ceramic 
composite.  Such  a  layer  may  develop  many  transverse  cracks  due  to  ther¬ 
mal  stresses  produced  by  manufacturing  cooldown,  fig.  2.  These  cracks  are 
roughly  orthogonal  in  fiber  directions  of  the  layers.  In  a  ceramic  fiber  com¬ 
posite  laminate,  for  example  a  SiC  matrix  reinforced  by  graphite  fibers,  the 
stiflhess  of  the  interphase  layer  is  of  the  order  of  the  stiffness  of  the  layer 
material.  For  large  crack  density  in  the  layer  the  shear  moduli  decrease  very 
significantly,  but  not  G„„.  Therefore  such  a  cracked  layer  can  be  considered 
as  an  interface  which  is  perfect  for  normal  contact  ,  [u„]  =  0  ,  but  is  im¬ 
perfect  for  shear.  In  that  event  the  term  in  G22,  equ.(17),  containing  D„, 
is  negligible. 
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CRACKED  INTERPHASE 

Figure  2.  Laminate  with  damaged  interphase 


3.3.  THERMAL  EXPANSION 

To  evaluate  the  effective  thermal  expansion  coefficient  (TEC)  of  a  cracked 
laminate  it  is  very  convenient  to  use  the  Levin  relation  as  has  been  done 
for  perfect  interlaminar  Interface  in  Hashin  (1988).  For  this  purpose  con¬ 
sider  any  elastic  composite  which  is  subjected  to  the  homogeneous  traction 
boundary  conditions  (7)  and  let  the  internal  stresses  due  this  loading  be 
a^(  as).  Denoting  the  local  TEC  a(*)  and  the  effective  TEC  o:*,  the  Levin 
relation  ,  Levin  (1967),  is  expressed  as 

[  cc.<r^dV=  a*.a^V  (24) 

Jv 

Levin’s  original  derivation  of  (24)  is  based  on  displacement  continuity,  but 
it  may  be  shown  that  it  remains  valid  for  interface  displacement  continu¬ 
ities  which  obey  the  relations  (3)  and  therefore  (24)  may  be  employed  in 
the  present  case  with  the  stresses  (9-10)  and  the  functions  ipm  to  give  an  ap¬ 
proximate  expression  for  the  TEC.  For  the  loading  AT^x  the  only  surviving 
component  of  in  (7)  is  =  o®  as  defined  by  (15).  Then  from  (24) 


rV  = 


dx  (25) 


L  is  the  length  of  the  laminate  of  unit  thickness  in  y  direction.  In¬ 
sertion  of  the  stresses  into  (25)  with  use  of  the  boundary  conditions  of  (16) 
yields 


oc. 


=  «L  + 


tri  1 

oq  1  +  A 


(a^  -  ay)  <  ^  > 


(26) 
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where  is  the  TEC  of  the  uncracked  laminate,  otx  and  ar  are  the  axial 
and  transverse  TEC  of  the  unidirectional  fiber  composite  and  <  ^  >  is  the 
average  of  the  random  variable  <j)^  which  is  defined  as 


3.4.  EFFECTIVE  SHEAR  MODULUS 

Let  the  cracked  laminate  shown  in  fig.  1  be  subjected  to  constant  shear 
membrane  load  which  defines  the  average  appUed  shear  stress 

T®  =  N^/2h  (29) 

In  this  case  the  laminate  is  in  a  state  of  antiplane  stress  with  respect  to  the 
y  axis.  Admissible  stresses  are  defined  as  in  Hashin  (1985)  by 

Then  a  variational  optimization  as  done  above  and  in  Hashin  (1985)  yields 
the  results 


'I’mi.O  —  Cosh{(l^)/Cosh{lipm) 

Ga 

-  1+  <  Tanh{pp)  >/Xp<p> 

In  the  case  of  equidistant  cracks  =  a  ,  Pm  —  aAii  and  all  of  the 
averages  in  all  of  the  expressions  above  reduce  to  simple  functions  of  p 
which  are  defined  by  removal  of  the  brackets. 

3.5.  STRESS  ANALYSIS  AND  CRACK  OPENING  DISPLACEMENTS 

The  optimal  functions  4>m  and  V’m  define  optimal  admissible  stresses  by 
the  relations  (9,10,30).  The  effective  properties  and  G^y  based  on  these 
are  strict  lower  bounds  which  also  agree  quite  well  with  experimental  data. 
The  status  of  the  stresses  associated  with  the  optimal  functions  is  less  clear 
which  is  a  typical  situation  for  any  variational  field  approximation.  It  is, 
however,  believed  that  these  stresses  are  of  qualitative  importance,  at  least. 
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and  the  numerical  results  obtained,  some  of  which  are  shown  below,  support 

this  belief.  x-  4.  *1, 

The  results  obtained  in  this  work  can  be  easily  used  to  estimate  the 

crack  opening  displacements  (COD)  when  the  cracks  are  equidistant.  It  is 
rigorously  true  that  the  stress  energy  17  of  a  cracked  elastic  body,  homoge¬ 
neous  or  non- homogeneous,  and  the  stress  energy  Uq  of  the  same  uncracked 
body,  under  same  load,  are  related  by 

+  (32) 

2  ^JSm 

where  is  the  surface  of  the  crack,  r°  is  the  traction  on  same  surface 
in  the  uncracked  body  and  [u]  is  the  COD.  In  the  case  of  simple  tension 
discussed  above  U  is  given  by  (15)  and  Uq  =  I2E^)V  for  the  uncracked 
laminate.  Also,  the  only  surviving  Tf  is  =  ui.  The  COD  is  now  estimated 
in  the  form  of  two  equal  and  oppositely  joined  second  order  parabolas.  Thus 

M  =  2tf(l-(|-)"] 

where  26  is  the  maximum  COD.  It  then  follows  easdy  that 

^Ai  =  ^°(1/^:  -  (33) 

2, 


4.  Results 

Illustrative  results  are  presented  for  a  [0“,  90°],  laminate  in  which  the  layers 
of  equal  thickness  are  T300/SiC  ceramic  unidirectional  composites  with 
fiber  volume  fraction  0.45.  The  relevant  properties  of  the  layer  material 

Ba  =  431.5  GPa  Er  =  113.6  GPa 

Ga  =  90.8  GPa  Gt=  39.3  GPa 

t/A  =  0.182  1^  =  0.446 

=  2.39  10-e(G°)  ay  =  5.49  lO-e(GO) 

It  is  assumed  that  in  between  the  layers  there  is  a  thin  oxidation  protec¬ 
tion  interphase  of  isotropic  P4G  material  with  thickness  0.02  of  the  layer 
thickness  and  with  properties 

Ei  =  380  GPa  Gi  =  159.7  GPa  i/,-  =  0.19 

Due  to  thermal  treatment  the  interphase  may  develop  many  through 
cracks.  As  has  been  explained  above,  this  creates  an  orthotropic  interphase 
which  may  be  considered  perfect  in  normal  z  direction  but  imperfect  in 


12 


shear.  Increase  of  crack  density  of  the  interphase  may  be  considered  as 
decrease  of  the  effective  shear  stiffness  of  the  interphase.  All  of  the  results 
given  below  are  in  the  form  of  plot  families  where  each  plot  is  associated 
with  a  shear  modulus  value  Gi/m  ,  m=l,5,  10,  20,  50,  100,  200.  Figs.  3—5 
show  plots  of  Young’s  modulus,  TEC  and  shear  modulus,  all  normalized 
with  respect  to  their  values  for  the  uncracked  laminate,  as  functions  of  crack 
density  for  the  case  of  equidistant  cracks,  expressed  by  the  parameter  p.  For 
large  values  of  p  the  properties  of  the  uncracked  laminate  are  attained  while 
for  small  values  of  p  the  properties  reduce  to  those  of  a  laminate  in  which  the 
90®  layer  has  vanishing  Ej*  and  Gj^  but  retains  it’s  value.  Such  stiffness 
loss  is  associated  with  the  concept  of  laminate  netting  analysis.  The  values 
of  the  properties  decrease  with  decreasing  interphase  shear  modulus.  Thus 
for  each  property  the  uppermost  plot  is  for  undamaged  interphase  which 
may  be  regarded  as  a  perfect  interface  while  the  lowest  plot  is  for  m=200.  It 
is  seen  that  the  effect  of  interface  imperfection,  i.e.  interphase  damage,  is  not 
very  significant  for  Young’s  modulus  and  TEC  but  is  very  significant  for  the 
shear  modulus.  It  is  also  seen  that  the  plots  for  normalized  Young’s  modulus 
and  TEC  are  very  similar  and  indeed  these  normalized  quantities  are  almost 
the  same  numerically.  It  should  be  realized  that  interphase  damage  has  no 
effect  on  an  uncracked  laminate  under  in-plane  loading  since  there  are  no 
interlaminar  stresses  in  this  case. 

Figs.  6-7  show  internal  stresses  as  functions  of  x  ,for  load  JV**,  when 
the  intercrack  distance  is  2a  =  5ti.  Fig.6  shows  the  in-plane  stresses  cf^x  iu 
the  cracked  layer  and  in  the  uncracked  layer  functions  of  x  for  the  case 
when  the  intercrack  distance  is  2a  =  5ti,  the  family  of  plots  being  defined 
by  the  sequence  Gi/m.  It  is  seen  that  the  tensile  stress  in  the  cracked 
90®  layer  decreases  with  increasing  interphase  damage  amd  therefore  the 
stress  in  the  uncracked  0®  layer  increases  with  damage.  Fig.  7  shows  similar 
plots  for  the  interlaminar  stresses  Cizi^th)  a-nd  crxz(x,ti).  These  stresses 
decrease  with  increasing  interphase  damage  and  the  effect  is  significant. 


S.  Effective  shear  modulus  versus  crack  spaang 


Normalized  In- Plane  Stresses 


Figure  6.  lu-plane  stresses  in  layers 
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•Normalized  Interface  Stresses 


5.  Conclusion 

Previously  developed  variational  analysis  of  cracked  laminates  has  been 
extended  to  the  case  of  imperfect  interlaminar  interface  by  use  of  a  gen¬ 
eralized  thermo-elastic  variational  principle  for  imperfect  interface.  In  the 
present  work  detailed  analysis  has  been  confined  to  a  cross-ply  laminate 
in  which  only  the  90®  layer  is  cracked.  But  there  is  no  difiiculty  to  carry 
out  similar  analysis  with  imperfect  interface  for  the  case  when  all  layers 
are  cracked,  on  the  basis  of  the  admissible  stress  system  which  has  been 
constructed  in  Hashin  (1987)  for  orthogonally  cracked  laminates.  It  is  also 
a  straightforward  matter  to  analyze  mechanical  and  therm^  stresses.  ^ 
Present  analysis  has  shown  that  for  a  cracked  laminate  with  interlaminar 
interface  which  is  imperfect  in  shear  only,  the  quantitative  effect  of  interface 
imperfection  on  effective  Young’s  modulus  and  TEC  is  not  drastic,  but  there 
is  significant  effect  on  the  effective  shear  modulus  and  internal  interlaminar 
stresses. 
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DIFFERENTIAL  SCHEME  FOR  EFFECTIVE  THERMOELASTIC 
PROPERTIES  OF  CRACKED  COMPOSITE  MATERIALS 


V.  Vinogradov  and  Z.  Hashin 


Differential  Scheme  for  Effective  Thermoelastic 
Properties  of  Cracked  Composite  Materials 

Introduction 

Analytical  determination  of  effective  properties  of  composite  materials  containing  ran¬ 
domly  oriented  microcracks  becomes  a  much  more  difficult  problem  than  the  problem  of 
small  crack  density  discussed  in  the  previous  chapter.  An  undefined  microgeometry  and 
interaction  between  constituents  make  the  simple  analytical  approach,  like  in  the  above- 
mentioned  case,  to  be  inapplicable.  Thus,  we  must  resort  to  an  approximate  scheme. 
In  this  chapter  the  differential  scheme  (DS)  approximation  will  be  applied  to  predict 
effective  thermoelastic  properties  of  composite  materials  containing  microcracks. 

The  DS  is  based  on  the  notion  of  incremental  construction  of  the  composite  material 
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by  gradual  addition  of  infinitesimal  amounts  of  inclusions.  According  to  this  idea,  for  a 
two- phase  composite  of  material  1  with  moduli  and  material  2  with  moduli  \  we 
start  with  homogeneous  matrix  material  1.  Then,  inclusions  of  phase  2  are  embedded 
in  the  matrix  in  such  a  way  that  the  inclusions  are  in  dilute  suspensions.  This  results 
in  effective  moduli  C*  which  may  be  obtained  by  the  well-known  solution  for  dilute 
concentration  of  inclusions.  The  next  stage  consists  of  imbedding  grains  of  phase  2, 
that  are  in  order  of  magnitude  larger  than  the  previous  ones  (so  that  the  fundamental 
assumption  that  they  "see”  an  effective  medium  with  moduli  C*  is  quite  accurate),  and 
so  on.  The  DS  hierarchy  of  composites  continues  until  phase  2  occupies  its  correct  volume 
fraction. 

We  shall  describe  a  procedure  whereby  a  JV-phase  composite  can  be  realized.  The 
first  of  the  two  procedures  considered  here  is  characterized  by  keeping  the  volume  of  the 
composite  fixed  at  during  the  entire  process,  called  the  Fixed  Volume  Process  (FVP),  and 
in  the  second,  the  volume  of  the  composite  is  allowed  to  increase,  the  Variable  Volume 
Process  (VVP).  A  schematic  representation  of  the  processes  is  shown  in  Fig.3.1.  It  will 
be  shown  that  both  of  the  processes  lead  to  the  same  canonical  equations  for  the  effective 
moduli. 

Let  us  begin  with  volume  Vq  of  material  0,  an  anisotropic,  homogeneous,  linear  elastic 
solid  with  elastic  moduli  tensor  which  we  refer  to  as  the  ’’backbone”  material. 

Grains  of  materials  1, 2, . . . ,  iV  are  imbedded  in  material  0  in  such  a  way  that  the  volume 
remains  fixed  at  Vo-  This  involves  removing  and  replacing  some  of  the  original  material 
(FVP).  Materials  1, 2, . . . ,  AT  are  also  anisotropic,  homogeneous,  linear  elastic  with  moduli 
k  —  1,2, . . . ,  AT.  The  composite  now  has  a  new  modulus  tensor  C*,  different  from 
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Figure  3.1:  The  iterative  construction  process  of  the  differential  scheme 
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The  construction  process  continues  by  removing  the  current  material  and  replacing 
it  with  grains  of  the  N  materials  until  the  phase  volume  fractions  Vk  reach  their  desired 
values 

Let  uo,  wi,  V2,...,VN  be  the  current  volume  fractions  of  the  materials  0, 1,2, 
respectively,  such  that 

Vo  +  +U3  +  . . .  +  Vat  =  1  (3.1) 

At  each  replcicement  the  removed  material  must  have  the  same  volume  fractions  of  the 
materials  0, 1, 2, . . . ,  AT  as  in  the  total  volume.  Thus,  at  each  stage  the  material  is  assumed 
to  be  homogeneous.  This  can  be  realized  if  the  replacement  grains  are  an  order  of 
magnitude  greater  in  size  than  those  at  the  previous  replax:ement.  In  addition,  the  grains 
must  be  dispersed  at  random  and  occupy  an  infinitesimal  volume  fraction. 

Let  the  process  be  parameterized  by  an  artificial  parameter  t,  which  is  analogous  to 
time.  The  origin  of  the  time  is  conveniently  chosen  so  that 

ufc(<  =  0)  =  0,  =  (3.2) 

For  each  t  the  material  is  homogeneous,  i.e.  any  unit  volume  of  the  composite  contains 
Vk{t)  of  material  k. 

The  volume  AF(t)  of  the  material  removed  at  any  time  t  is  replaced  with  volumes 
AVk{t)  such  that 


AVk{t)  +  AV2{t)  +  . . .  +  AyA^(t)  =  AVk{t)  =  AV(t) 


k=i 


(3.3) 
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After  the  replacement  the  volume  of  material  k  is 


N 


V  [t;fc(t)  +  At;jk(t)]  =  V Vk{t)  -  Vk{t)^AVi{t)  +  AVk{t) 


i=l 


(3.4) 


where  V  =  is  the  volume  of  the  composite,  the  first  term  in  the  right  part  of  the 
equality  is  the  volume  of  phase  k  before  the  replacement,  the  second  term  defines  the 
removed  volume  and  the  third  one  represents  the  added  volume  of  phase  k.  This  relates 
the  increments  in  volume  fractions  Vk{t)  to  the  corresponding  increments  in  Vk{t).  When 
AV{t)  is  small,  these  equations  assume  the  form 


N 


t=l 


(3.5) 


where  /(<)  =  -i-ZCO*  This  relation  between  Vk{i)^  ^k{t)  and  Vk{t)  can  be  inverted  to  give 

at 


Vk{t)  as; 


(3.6) 


where  v{t)  =  ^  ei(0>  6(f)  =  S  Vk{t). 

t=l  t-1 

For  the  VVP  process  the  same  relation  may  be  obtained  taking  into  account  that  V 

is  allowed  to  increase.  Begin  at  t  =  0  with  an  initial  volume  Vq.  At  each  step  of  the 

homogenization,  incremental  volumes  AVfc(t)  of  materials  1,2, ...,A  are  added  to  the 

N 

current  total  volume  y(t)  =  14  +  ]C  Vk{t)-  The  volume  fraction  of  material  k  is  defined 


as 


Vk{t)  = 


Vkjt) 

V{t) 


(3.7) 


Now  we  are  retaining  all  the  added  material  within  the  mixture.  The  derivative  of  (3.7) 
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with  respect  to  time  t  ha^  the  same  form  as  for  the  FVP: 


.  Vt(t)V(t)-v^(t)V(t)  v,{t) 


Consequently  the  expressions  for  Vk{t)/V{t)  are  the  same  as  (3.6) 


W) 

V{t) 


=  Vk{t)  +  Vk{t) 


^(0 

1  -  v{t) 


with  the  only  difference  that  the  constant  V  in  (3.6)  is  replaced  here  with  V{t),  which  is 
a  function  of  t. 

When  the  small  amount  of  inclusions  AF (<)  is  added  to  the  homogeneous  material 
with  stiffness  C*(t)  and  compliance  S*(t)  the  incremental  change  in  the  moduli  is  given 

by 

AC-(<)  =  x;  (cW-C-(<))  AW(()^(i) 

AS-(«)  =  i;(sl‘)-S-(())BW(t)^(() 

fc=l 

where  A(*^(t)  and  are  forth  rank  concentration  tensors  evaluated  for  dilute  con¬ 

centration  of  the  inclusions  in  the  matrix  with  moduli  C*(t),  and  defined  by  the  linearity 
relations 

Here  and  ^  are  average  strain  and  average  stress  over  phase  k,  the  first  when  e  is 
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prescribed  and  the  second  when  a  is  prescribed.  In  the  limit  as  t  >  0  we  have 


C-(()  =  E  (C<‘)-C'(i))  A(‘)(()h(t) 

k=l  / 

S-(<)  =  l;  (S“l-S'(t)) 

A:=l  ^ 


(3.8) 


where  ^(t)  are  defined  by  (3.6).  Because  the  relation  ^(<)  is  common  for  both  the 
FVP  and  the  WP,  differential  equations  (3.8)  are  said  to  be  canonical  equations  of  the 
DS,  independent  of  the  specific  material  replacement  process.  Integration  of  (3.8)  with 
use  of  initial  conditions  C*{t  =  0)  =  and  S*(t  =  0)  =  gives  the  effective  elastic 
moduli  of  multi-phase  composite  material. 

It  must  be  noted,  that  functions  Vk{t)  describe  a  path  in  the  (vi,V2,...,UAr)  space 
as  t  varies,  and  t  may  be  considered  as  an  arc  length  parameter.  The  path  begins  at 
the  origin,  it  is  restricted  by  the  plane  ui  -1-  W2  "b  •  •  •  "b  1  and  terminates  at  the 


point  (v(,V2,..  •  jujv),  which  is  associated  with  volume  the  fractions  of  the  phases  in  the 
composite  under  consideration.  An  example  of  such  a  path  for  a  three-phase  composite 
is  shown  in  Fig.3.2.  It  is  apparent  that  any  arbitrary  functions  ujt(f),  which  lead  to  the 
same  path  in  the  (ui,U2)-  •  •  jViv)  space,  will  give  equal  results  for  the  effective  moduli. 
Additional  relations  between  the  variables  Vk{t)  and  Vk{t)  for  the  construction  procedures 
described  above  are  given  in  NORRIS  [20]  and  NORRIS  et  al.  [21].  An  explanation  of  the 
distinction  between  paths  will  be  given  in  section  3.3. 
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0  v|  1 

Figure  3.2:  A  homogenization  path  in  the  (v  j  .v  j  )  plane. 
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3.2  Effective  Thermal  Expansion  Coeflacients 

3.2.1  Development  of  Effective  Thermal  Properties 

Consider  a  heterogeneous  body  of  volume  V  subjected  to  temperature  change  <p  (Fig.3.3). 
The  body  consists  of  an  arbitrary  anisotropic  matrix  0  with  moduli  and  thermal 
expansion  coefficients  and  inclusions  with  moduli  =  1,2, .. .  ,-A^)  and  thermal 

expansion  coefficients  oc  .  At  each  point  of  the  body  the  stresses  and  strains  are  related 

by 

e  (x)  =  S(x)<r  (x)  +  a  (x)(^ 

cr(x)  =  C(x)e(x)  +  r(x)si? 

where  r(x)  =  -C(x)a  (x)  and  S(x)  =  C“^(x)  is  the  compliance  tensor.  For  given  bound¬ 
ary  conditions  (T(5)  =  0;  (p)  the  average  stresses  and  strains  are  defined  by 

J  / ,(x)dV  =  i  /  e(x)dV+^  J  4K)dV=  a  V 

^  V  AV  V-AV  (3 10) 

j cr{:<.)dV  =  ^  J  <r(x)(iF  +  —  j  o-(x)dV  =  0 

^  V  iv  V-AV 


Figure  3.3:  A  composite  body  under  temperature  change  (p  {T{S)  =  0) . 
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where  <x*  are  the  effective  thermal  expansion  coefficients  of  the  composite,  AV  defines 
the  volume  of  inclusions  and  V  —  AV  defines  the  volume  of  the  matrix.  After  substitution 
of  (3.9a),  the  last  equality  in  (3.10a)  may  be  rewritten  in  the  form 

=  J^i^)dV+  j  e{x)dV 

AV  V-AV 

=  J  [S(x)o-(x) +  a(x)v?]dV'+  J  [S^°V(x)  + 

AV  V-AV 

and  by  using  the  fact  that  the  average  stress  is  equal  to  zero  (3.10b)  we  obtain  the 
following  result: 

c.*ipV  =  cc^\V  +  f  [S(x)-S(°)]<r(x)dV^+  J  [a(x)-a(”)]<^dy  (3.11) 

AV  AV 

Because  the  compliance  tensor  S(x)  is  a  piecewise  constant  function,  the  integration  in 
this  equation  may  be  replaced  by  the  sum: 

+  l;  [sW-SW]  i  j  <.{x)dV  +  I;  (3.12) 

k=l  k=l 

Now  recalling  the  definition  of  average  stress,  this  equation  may  be  rewritten  in  the  form 

=  <.<»v  +  E  [S(‘)-S(“)]  J ^  +  E  (3.13) 

it=l  ^  ^  k=l  ^ 

or  eliminating  the  temperature  (p  we  obtain  the  exact  result 

A«  =  E  [sW-S*"*]  +  E  ^  (3.14) 

fc=l  ^  fc=l  ^ 
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where  Aa  =  a  *  —  a  and  is  a  second  rank  tensor  defined  as  the  average  stress  in 
phase  k  due  to  an  applied  unit  temperature  change: 

=  (3-15) 

Now  return  to  the  DS  described  in  the  previous  section  and  consider  the  DS  con¬ 
struction  process  at  any  instant  of  time  t.  An  addition  of  inclusions  of  volume  AV  to 
the  homogeneous  material  with  overall  moduli  C*  and  TEC  a*  will  lead  to  the  change 
of  overall  TEC  Aa*,  which  according  to  (3.14)  assumes  the  form 

=  E  [sW-S'W]  ^(()  +  E  (3-lS) 

In  the  framework  of  the  DS,  at  each  step  of  the  construction  process  only  infinitesimal 
amount  of  inclusions  is  embedded  in  the  homogeneous  material  obtained  at  the  previous 
step.  Thus,  equation  (3.16)  may  be  parameterized  by  the  artificial  time  t 

«•(*)  =  E  [s<‘'-s-(()]  ^(«)  +  E  vW  (317) 

Jt=l  ^  A;=l  '' 

and  tensor  b(*'^(i)  may  be  obtained  by  solving  the  problem  where  dilute  concentration 
of  inclusions  of  phase  k  is  embedded  into  the  matrix  with  effective  properties  C*{t)  and 

The  differential  equations  (3.17)  may  be  simplified  by  resorting  to  the  so-called  decom¬ 
position  scheme,  proposed  by  BenvenistE  &  DVORAK  [23].  According  to  their  result, 
the  concentration  tensor  b^^^  can  be  uniquely  determined  in  terms  of  the  fourth  rank 
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concentration  tensor  which  relates  the  stresses  in  phzise  k  to  the  external  applied 
homogeneous  traction  Ti{S)  =  ,  and  is  defined  as 


Their  result  can  be  rewritten  for  our  case  in  the  form 


where  I  is  the  fourth-order  unit  tensor.  Substitution  of  this  relation  into  (3.17)  leads  us 
to  following  differential  equations  for  components  of  the  effective  TEC  tensor 


«*(t)=  f:  -  S*(<))  bW(0  -  S*(t))  ' 

fc=l 


(3.19) 


when  ^(0  defined  by 


14 

V 


(t)  =  Vk{t)  +  Vk{t) 


1  -  v{t) 


Equations  (3.19)  and  (3.8b)  and  the  additional  initial  conditions  a*{t  =  0)  =  give  a 
set  of  coupled  differential  equations  for  estimation  of  the  effective  thermoelastic  moduli 
of  multi-phase  composites. 

Now  assume  alternative  boundary  conditions,  when  homogeneous  displacements  u(S')  = 
e^x  =0  are  applied  to  the  external  surface  of  the  body  (see  Fig.3.4).  In  this  case  the  ex¬ 
pressions  for  e  and  a  becomes 
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Figure  3.4:  A  composite  body  under  temperature  change  («(5)  =  0) . 

«-=i/.(x)<iV'  =  0 

(3.20) 

^  j .,(-x)dV  , 

V 

Analogous  consideration  of  the  problem  results  in  the  following  set  of  coupled  differential 
equations  for  components  of  tensor  T“ 

E  W)  (r('=)-r*(<))  y{t)  (3.21) 

k=zl 

which  together  with  (3.8a)  and  the  initial  conditions  r’(t  =  0)  =  may  be  used  to 
predict  the  overall  thermal  stress  tensor  F*. 

3.2.2  Proof  of  Consistency 

It  is  required  to  prove  that  ot*  and  F*  cls  predicted  by  (3.19)  and  (3.21)  are  mutually 
reciprocal  in  the  sense  that 

F*=  -C*a*  (3.22) 
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If  the  effective  themoelastic  moduli  obey  this  relation  at  the  endpoint  of  the  DS  homog¬ 
enization  process,  it  must  also  be  satisfied  at  each  instant  t 


r-(()=  -C>(t)a‘(() 


Let  us  take  the  derivative  of  r*(t)  with  respect  to  time  t 

^r-(()=  (c-(t)o-(f))  =  -^c*(i)o.-(t)  -  c-(i) 

Substituting  expressions  (3.19)  and  (3.8a)  for  <i*(<)  and  C*{t)  respectively  into 
obtain 


r-(*)  =  -i;(cW-c-(t))AW(()<.*(Ov(*) 

k^l  ^ 

(S^''^-S*(t))  B^*)(t)  (S('=)-S*(f))”'  (a^'')-a*(t)) 

k=i 

which  after  simplification  of  the  right  side  may  be  rewritten  in  the  form 


f*(t) 


a('')(<) 

fc=l 


X  <«*(()-S*(<)  {s<‘’>-S‘(t)) 


where  the  relation  between  concentration  factor  tensors  B(^)(t)  and  A^*')(t), 
dilute  concentration  of  inclusions  as 


BW(t)  =  CWA(*^>(t)  S*(t) , 


(3.23) 

(3.24) 

(3.24)  we 

(3.25) 
j{t) 

(3.26) 

given  for 

(3.27) 
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was  used.  Recalling  the  equality  C*{t)  S*{t)  =  I  proved  for  the  DS  by  Hashin  [14]  and 


the  identity  that  results  from  it: 


s*(t)  (sW-s*(t))  '  =  -  (cW-c*(t)) 


differential  equations  (3.26)  reduce  to 


N 


Vk 


k-l  ^ 


which  faithfully  copies  (3.21)  obtained  by  direct  treatment  of  (3.20). 

This  proves  that  (3.23)  is  satisfied  and  completes  the  self-consistency  of  the  DS. 


3.2.3  Levin’s  Formula 


Consider  a  two-phase  material  which  consists  of  the  matrix  with  compliance  and 
thermal  expansion  coefficients  and  inclusions  with  compliance  and  TEC 
In  this  case  equations  (3.8b)  and  (3.19)  reduce  to 


s-(()  =  (s<')-s-(())bI'>(«)yz^ 

ci-(i)  =  (s<»-s'(())b<''(«)  s'w) 


Eliminating  the  common  multipliers  from  these  equations  we  obtain 


a*(f)  =  S*{t)  (S(’)  -  S*(t))  ^  («(’)-«*(t)) 


(3.28) 
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Now  define 


x(f)  =  q:*(<) 

Y(f)  = 


Accordingly,  the  derivatives  of  these  tensors  with  respect  to  t  are 


x(<)  = 


Y{t)  =  -S\t) 


and  equations  (3.28)  assume  the  form 


x(<)  =  Y(f)Y(()-'x(<) 


(3.29) 


The  solution  of  these  coupled  differential  equations  may  be  written  as 


x(t)  =  Y(t)  Y(f)  =  Y (t)  z{t) 


where  z{t)  =  Y{t)  ^x(<).  Substitution  of  this  into  (3.29)  leads  to 


Y{t)  z{t)  +  Y(t)  i{t)  =  Y(t)  Y(t)-'Y(0  z{t) 


Noting  that  the  expression  in  the  right  side  and  the  first  term  in  the  left  one  are  equal, 
we  obtain 


Y(^)  z(t)  =  0 
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which  results  in  the  following  simple  form  for  z(t) 


=  Cij=const. 

Therefore 

x(t)  =  Y(t)c  (3.30) 

and  recalling  the  definitions  of  x(t)  and  Y(t)  we  can  rewrite  (3.30)  in  the  form 

«*(<)  =  «  -  S*(<))  c  (3.31) 

where  Cij  may  be  foimd  from  the  initial  conditions  at  f  =  0: 

cc*{t  =  0)  =  «(®) 

S*(t  =  o)  =  s(°) 

It  follows  that 

C  =  -(SW-S<'))'‘(a  <’)-«»>) 

which  together  with  (3.31)  yields  the  relation  between  a*(t)  and  S*(t)  at  any  instant  t 

«■(()  =  «<')  +  (s<')  -  S-(i))  (s«>  -  (al')-aP)) 

and  consequently  at  the  endpoint  of  the  homogenization  process 

a*  =  «(l)  +  -  S*)  (S(2)  -  (3.32) 
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This  result  is  the  well-known  Levin’s  formula  (see  LEVIN  [24]).  It  is  interesting  to 
note,  that  the  result  of  the  DS  approximation  satisfies  Levin’s  formula,  which  is  precise  for 
any  two-phase  statistically  homogeneous  composite  material  with  completely  arbitrary 
phase  geometry  and  elastic  symmetry. 


3.3  Path  dependence  of  the  homogenization  process 

The  coupled  systems  of  differential  equations  (3.8a),(3.8b),(3.19)  and  (3.21)  depend  upon 
the  path  we  choose  to  consider  in  the  (ui,U2,  •  •  •  j  vn)  space.  In  fact,  two  different  paths 
from  the  origin  to  the  same  endpoint  (v{^V2, . .  • ,  vjv)  will,  in  general,  not  give  the  same 
result  at  the  endpoint.  In  this  section  an  effort  will  be  made  to  explain  the  distinction 
between  paths. 

The  key  to  the  present  consideration  is  the  idea  of  realization,  introduced  by  MiLTON 
[25].  An  approximate  scheme  is  called  realizable  if  we  can  describe  a  construction  process 
to  make  a  composite  material  with  effective  moduli  predicted  by  the  scheme.  Thus,  the 
Hashin-Shtrikman  bounds  for  the  bulk  modulus  of  a  two  phase  randomly  disordered  solid 
can  be  realized  through  the  well  known  packed-spheres  geometry  (e.g.  Hashin  [26]). 

Let  us  list  the  features  of  the  inner  structure  of  composite  materials,  obtained  by  the 
DS  construction  processes  described  above. 

1)  The  essential  characteristic  is  the  diversity  of  inclusion  sizes. 

2)  Inclusions  of  each  size  are  in  dilute  concentration, 

3)  They  are  dispersed  at  random  and  sufficiently  separated. 
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The  first  feature  suggested  by  ROSCOE  [27]  arises  from  the  fact  that  at  each  construc¬ 
tion  step  the  composite  ha^  to  look  homogeneous  to  the  inclusions  being  inserted,  which 
is  possible  only  if  the  inclusions  are  increasing  in  size  at  each  succeeding  step  of  the  con¬ 
struction  process.  More  specifically,  the  size  of  inserted  inclusions  must  be  comparable 
to  the  representative  volume  element  (RVE)  size. 

From  this  point  of  view,  different  paths  in  the  (wi,U2>-  •  •  j^jv)  space  may  be  associ¬ 
ated  with  different  microstructures  of  a  composite,  or  in  other  words,  with  different  size 
distributions  of  inclusions. 

Consider  a  composite,  which  satisfies  the  above-enumerated  requirements  for  material 
microstructure.  A  composite  of  this  kind  contains  inclusions  of  different  sizes  Zi,  which 
may  be  ordered  as  the  increasing  sequence 

Z-l,  Z2,  ^3,  •  •  •  ,  ■  •  .  ,  (Af  OO) 


SO  that 


Zi  Z2<.  Zz  Zr  ZM 


Each  group  of  inclusions  characterized  by  a  common  size  Zj.  may  contain  a  number  of 
phases.  The  volume  fraction  of  phase  k  in  group  r  is  denoted  Au{(r).  In  this  case  the 
following  relations  are  satisfied 


M 


5^Auf(r)  =  vi 

r=l 

Y^Avlir)  =  Av^{r) 


k=i 
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where  Av^  (r)  is  the  volume  fraction  occupied  by  group  r. 

According  to  this,  we  can  build  the  construction  process  starting  with  inclusions 
of  size  zi,  afterwards  embedding  inclusions  of  size  Z2  and  so  on,  up  to  zm,  where  the 
subscript  index  denotes  now  the  step  of  the  construction  process. 

First  consider  the  WP,  according  to  which  all  the  volume  inserted  to  the  composite 
remains  till  the  end  of  the  construction  process  and  creates  the  final  microstructure. 
Thus,  using  formula  (3.7)  we  obtain  the  points  Ufc(r)  on  the  homogenization  path 

Vkir)  =  — -  (3.33) 

t=i 

An  example  of  a  homogenization  path  built  for  a  given  3-phase  microstructure  is 
shown  in  Fig.3.6.  Size  distribution  of  the  inclusions’  volume  fractions  of  is  shown  in 
Fig.3.5a  for  phase  1  and  in  Fig.3.5b  for  phase  2.  The  homogenization  path  obtained  by 
(3.33)  is  represented  in  Fig.3.5c. 

The  same  equation  (3.33)  may  be  obtained  from  consideration  of  the  FVP,  when  the 
volume  of  phase  k  added  at  step  r,  Aufc(7')j  decreases  repeatedly  during  the  rest  of  the 
construction  process.  For  any  phase  k  the  following  relationships  are  satisfied. 

1.  The  volume  added  at  the  last  step  M  is  not  changed: 

Avk{M)  =  Avl{M) 

2.  The  volume  added  at  step  (M  —  1)  decreases  because  of  the  removing  of  Avl(M) 
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r,  number  of  particle  size 
or  step  of  constr.  process 


V,  =  0.4 


r,  number  of  particle  size 
or  step  of  constr.  process 


at  step  M: 

Avl{M  -  1)  =  Avk  {M  -  1)  (1  -  Av{M)) 


3.  The  same  as  (2)  for  any  step  r: 


M 

Au{(r)  =  Aufc(r)  ]][  (1  -  At;(t))  (3.34) 

{=r-f  1 


On  the  other  hand,  the  volume  fraction  An*  (m)  of  the  inclusions  characterized  by  size 
that  are  added  at  any  step  m  (m  <  r)  will  decrease  towards  the  step  r,  and  will  be 
equal  to 

Awjt(m)  f[  (l-Av(z)) 

Combining  these  expressions  for  m  =  1  to  r  we  obtain  the  volume  fraction  of  phase  k  at 
step  r  of  the  construction  process 


=  XI  ("*)  n  (1  “  Au(z))  (3.35) 

m=l  z=:m+l 


Substituting  (3.34)  into  (3.35)  we  obtain 


(”i) 

n  (1-At;(i))  n  (1-At;(i)) 

t=m*fl  J  t=r+l 


Using  again  the  recursion  formula  (3.34),  the  product  in  the  denominator  of  this  equation 
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can  be  obtained  in  terms  of  the  group  volume  fractions  at  the  endpoint 

M  Mr 

II  (1-Av(t))  =  l-  2  (i)  =  ^ (i) 

t=r+l  t=r+l  *=1 

and  eventually 

vk(r)  =  - 

Vo +  (0 

which  is  precisely  equivalent  to  formula  (3.33)  for  the  VVP. 

This  deduced  relation  proves  that  the  homogenization  path  in  the  (ui,U2, . . . 
space  is  uniquely  determined  by  the  (st'^e  of  inclusions)— {volume  fraction  of  phase) 
dependences,  which  describe  the  microstructure  of  the  composite  imder  consideration. 
The  path  is  independent  of  choice  between  the  construction  processes  discussed  above, 
by  means  of  which  the  DS  is  interpreted. 

It  must  be  noted,  that  in  order  that  the  DS  equations  give  a  reasonable  approximation 
for  elFective  moduli  of  a  composite,  its  inner  structure  must  comply  with  the  above-hsted 
requirements  for  the  class  of  microgeometries  obtained  by  the  DS.  If  this  is  not  the  case, 
or  if  the  microstructure  is  not  defined  and  only  volume  fractions  of  phases  are  known, 
some  assumptions  of  the  path  shape  must  be  made  to  estimate  the  effective  moduli  of 
the  composite  material.  For  example,  in  the  case  when  all  inclusions  have  the  same  size, 
the  straight  line  from  the  origin  to  the  endpoint  (v{,v(,. . .  ,v^)  may  be  chosen  as  a 
homogenization  path. 

In  addition,  we  shall  emphasize  that  details  of  the  material  microgeometry  are  im¬ 
possible  to  obtain  from  a  given  path.  If  the  unique  microstructure-path  correspondence 
exists,  the  path-microstructure  dependence  does  not.  As  an  example,  the  straight  line 
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from  the  origin  to  the  endpoint  corresponds  to  an  infinite  number  of 

size- volume  distributions,  for  which  the  phase  volume  fractions  of  inclusions  of  each  size 
are  in  the  same  proportion  as  the  volume  fractions  of  phases  in  the  entire  composite. 

3.4  DS  for  Cracked  Composite  Materials 

The  purpose  of  this  section  is  to  apply  the  DS  approximation  to  estimate  the  effective 
thermoelastic  properties  of  composite  materials  containing  microcracks.  It  is  possible 
to  assign  crax:ks  to  a  special  class  of  inclusions,  and  thus  all  the  assumptions  about  the 
construction  procedure  can  be  extended  to  the  present  case.  At  eax;h  step  of  the  process 
infinitesimal  amounts  of  inclusions  and  cracks,  larger  in  size  than  at  the  previous  stages, 
are  added  to  a  homogeneous  material.  They  are  dispersed  at  random  and  well  separated 
in  such  a  way  that  the  material  remains  statistically  homogeneous  at  each  instant  t. 
Special  features  of  a  crack  as  an  inclusion  are  that  its  stiffness  vanishes  =  0,  and  it 
occupies  an  infinitesimal  volume  fraction  Vc  0. 

Let  us  rewrite  equations  (3.8a)  and  (3.19)  in  the  form 

=  E  (C<‘>-C-(<))  AW(i)  (3-37) 

JV+l  ^  .  _ 'y 

-  S*(t))  C('')A('=)(f)  S*(t)  -  S*(f)) 

fc=i 

X  (<.<‘>-«*(<))  (*.(()  +  (3-38) 

where  relation  (3.27)  was  used  and  phase  IV  +  1  represents  the  cracks.  Substituting 
=  0  into  these  differential  equations  we  obtain  that  the  N  +  1  term  in 
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(3.38)  vanishes,  and  (3.37)  reduces  to 


-C-(()A'(e)(i.(<) +  ».(«)  (3.39) 

N 

with  V  =  ^  Ufc  +  Uc.  Now  multiply  both  side  of  (3.39)  by  Vc  and  perform  the  limit 

k=l 

Vc(t)  — ♦  0  in  such  a  way  that  the  voids  flatten  out  into  cracks.  Then 
^  ^  ta  ^  [.,«)  E  (C<‘>-C-(())  aW«)  («(«)  + 

-c-WPW4^(”«w+''=Wr^) 

N  N 

where  t;  =  v*,  +  t>c  — X/  forth  rank  tensor  P(t)  is  defined  as 

k=l  k=l 

and  depends  on  moduli  C*(t)  of  the  material  around  the  cracks  added  at  time  t  and  the 
cracks’  geometry. 

As  an  example  we  specialize  (3.40)  to  the  important  case  when  the  cracked  particulate 
composite  material  and  its  constituents  are  statistically  isotropic,  Fig.3.6.In  this  event 
all  tensors  appearing  in  (3.40)  must  be  isotropic  and  we  represent  them  in  the  convenient 


Figure  3.6:  Schematic  diagram  for  randomly  distributed  cracks  and  inclusions 


form,  Hill  [28] 


C  =3  K  r+2  G  r=  (3/1',  2G) , 

A={Ar,A2), 

P=(A,P2), 

=  i  . 

iT=  r,  ri"  =  F,  rr  =  rr  =  o, 

=  ^kl  5  —  0 


(3.41) 


where  K  and  G  are  the  bulk  and  shear  moduli  respectively,  and  Ai,A2,Pi,P2  are  con¬ 
stants.  Thus  (3.40)  assumes  the  form 


Urn  G-  -  g  (G<‘I  -  G-)  (v,  +  »» j-^)  I  =  -G-Pi  Urn  (i,.  + 


(3.42) 


where  the  dependence  on  t  is  omitted  for  convenience. 


We  restrict  our  consideration  to  the  case  of  randomly  distributed  penny-shaped  cracks 
when  they  are  obtained  by  flattening  ellipsoidal  voids  which  have  two  equal  semiaxes  of 
length  a  and  the  third  one  pa,  where  p  =  const  can  be  made  indefinitely  small.  Then 
the  volume  fraction  Vc  is 


irp'£,a^ 

V 


=  rpa 


(3.43) 


where  a  without  indices  is  the  crack  density  parameter  (CDP)  introduced  in  Chapter  2. 
Substituting  this  expression  for  Vc  into  (3.42)  we  obtain 


JT-  =  g  (a-W  -  K-)  (4  +  (i  + 

<>  =  E  -  «•)  4‘>  (*.  +  -  <?•  t  (« + 

k=zl 


(3.44) 


Now  we  refer  to  the  results  of  Chapter  1  for  the  effective  bulk  and  shear  moduli  for  small 
density  of  open  microcrax:ks^(see  Table  2.1).  Assuming  that  a  is  small,  we  can  rewrite 
the  result  as 


K*-K  ^ 
G*-G  = 


-K 

-G 


9  (1  -  2u) 

32  (5  -v){l- 
45  (2  -  v) 


a 


which  for  any  step  of  the  DS  procedure  implies  the  form 


dK*  = 


-K* 


16  (1  - 1/*2) 
9  (1  -  2v*) 


da 


^The  load-induced  anisotropy  discussed  in  the  previous  Chapter  can  not  be  considered  in  the  frame¬ 
work  of  the  DS.  The  reason  for  this  is  that  if  the  applied  loads  change  the  material  symmetry  after 
the  first  step  of  cracks  addition,  the  successive  randomly  oriented  cracks  are  embedded  into  anisotropic 
material.  In  order  to  evaluate  tensor  P  we  have  to  know  the  crack  opening  displacement  [u]  (see  Chapter 
2)  for  a  crack  tilted  to  the  material  axes  of  symmetry,  which  is  an  extremly  complicated  problem. 
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dG* 


-G 


45  (2  - 1/*) 


Recalling  the  definition  of  tensor  P  it  may  be  concluded  that 


„  16  (1  - 

Pi  =  -X-,. - ::: — =  «  <3: 

9  (1  -  2u*) 

- - - - - (X  Z=z  TJ  a 

A  r'  /C\  _  .ifc\  t 


45 

with  I/*  = 


(2-1/*) 
ZK*  -  2G* 
2  {ZK*  +  (?•) 


(3.45) 


Therefore 


^  (iirw  -  K-)  4**  {vi + -  k‘k'  (a  + 

G-  =  f  (GW  -  G-)  4*’  (i.  +  -  GW  (a  +  a^) 


(3.46) 


where  k*  and  rj*  are  defined  by  (3.45). 

The  differential  equations  (3.38)  will  reduce  to 


=  E  K>  -  «•)  (i.  +  i^t;) 


(3.47) 


where  the  influence  of  the  cracks  enters  implicitly  through  K*. 

Again,  the  question  about  the  correspondence  between  composite  microstructure  and 
an  influence  path  is  raised  here.  The  relation  (3.33)  remains  unchanged  for  inclusions, 
and  for  cracks  it  assumes  the  form 


53Aq:^(0 

a(r)  =  - 

^0 
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where  Aa^ (i)  is  the  CDP  of  the  cracks  of  size  i  in  the  finite  material. 

As  an  example  consider  three  different  paths  to  the  same  endpoint  for  the  DS 

construction  process  by  which  a  three-phase  material  can  be  obtained  (Fig.3.7).  The  first 
one  is  the  straight  line  from  the  origin  to  the  endpoint  A  .  For  each  size  z,  the 

volume  fraction  of  the  particles  and  the  density  of  the  cracks  are  in  .the  same  proportion 
as  in  the  composite. 

According  to  the  next  path  OCA,  the  smallest  cracks  are  in  order  of  magnitude  larger 
than  largest  particles.  Thus,  the  construction  procedure  is  built  in  such  a  way  that  first 
the  inclusions  and  then  the  cracks  are  added  to  the  material.  The  first  stage  corresponds 
to  the  line  OC  along  the  axis  Ui,  and  the  line  CA  represents  the  insertion  of  the  cracks. 
Note,  that  CA  is  a  horizontal  line  because  the  embedding  of  the  cracks  does  not  change 
the  volume  of  the  composite,  and  hence  the  volume  fraction  of  the  inclusions. 

The  third  path  OB  A  is  opposite  to  the  previous  one.  We  start  with  small  cracks 

and  increase  the  crack  density  to  the  value - j  (point  B)  in  such  a  way  that  it  will 

\-v{ 

later  decrease  to  the  desired  value  oc^  (point  A)  during  removing  of  the  cracks  and  their 
replacement  by  the  larger  particles  (for  the  FVP),  or  during  the  material  volume  V 
increase  with  the  insertion  of  the  inclusions  (for  the  VVP). 

Elastic  moduli  K*  and  G*  and  thermal  expansion  coefficient  a*  of  the  cracked  com¬ 
posite  material,  containing  epoxy  matrix  and  40%  of  glass  spheres  as  inclusions,  are  cal¬ 
culated  for  the  described  paths  by  (3.46)-(3.47)  and  shown  in  Fig.3.8-3.10.  The  phase 
properties  are 
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E 

1/ 

a 

matrix 

3.45  GPa 

0.35 

42.  10"® 

particles 

72.40  GPa 

0.2 

5.  10-«°C-^ 

In  all  the  cases  of  Fig.3.8-3.10  we  show  the  effect  of  creiek  density  on  the  effective  prop¬ 
erties.  Both  K*  and  G*  vanish  asymptotically  with  increasing  crack  density.  Note,  that 
the  event  when  the  size  of  the  cracks  is  larger  than  the  particle  size  corresponds  to  the 
upper  curves,  and  the  opposite  case  to  the  lower  ones.  If  the  difference  in  the  results  for 
K*  due  to  the  influence  path  choice  is  insignificant,  the  shear  modulus  G*  depends  more 
heavily  on  the  path. 

The  dramatic  effect  of  the  path  on  the  effective  thermal  expansion  coefficient  a*  can  be 
seen  from  Fig.3.10.  If  small  cracks  and  cracks  that  are  comparable  to  particles  drastically 
decrease  the  value  of  a*,  the  cracks  which  are  significantly  larger  than  the  particles  have 
no  influence  on  the  effective  thermal  expansion  (the  upper  curve  in  Fig3.10).  This  can 
be  explained  by  that  these  cracks  ’’see”  an  effective  isotropic  material  which  consists  of 
the  matrix  and  the  inclusions  added  previously,  and  it  is  well-known  that  homogeneous 
heating  of  a  homogeneous  body  containing  cracks  does  not  produce  thermal  stresses  and 
does  not  change  the  thermal  expansion  of  the  body. 
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Figure  3.7:  Three  homogenisation  paths  in  the  (v,  CDP)  plane  to  the  same  endpoint  A. 
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Epoxy  matrix 
Glass  particles  40% 


Figure  3.8:  Variation  of  bulk  modulus  (GPa)  with  crack  density 


CDP 


Figure  3.9:  Variation  of  shear  modulus  (GPa)  with  crack  density. 
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MACRO-RESIDUAL  STRAINS  DUE  TO  CYCLIC  LOADING  OF 

COMPOSITES 


Z.  Hashin  and  B.W.  Rosen 


Macro-Residual  Strains  Due  to  Cyclic  Loading  of 
Composites 


INTRODUCTION 

It  is  well  known  that  when  a  metal  specimen  is  loaded  into  the  plastic  regime  and  then  unlo¬ 
aded  permanent  strain  which  is  known  as  residual  strain  will  remain.  This  phenomenon  is  expla¬ 
ined  by  internal  microscopic  changes  which  take  place  within  the  slip  systems  of  the  crystals  of 

which  the  metal  is  composed.  ,  -  r-  ,  *  , 

Similar  phenomena  are  observed  for  ceramic  matrix  composites  as  is  illustrated  in  fig.l  taken 
from  ref  [1]  Such  composites  consist  of  elastic-brittle  constituents  which  by  themselves  do  not 
exhibit  any  residual  strains  under  load  cycles.  Yet  when  such  composites  are  subjected  to  load 
cycles  residual  strains  on  the  macroscale  are  observed.  Characteristically,  the  residual  strains  grow 
with  the  number  of  cycles  until  they  attain  a  limit  value  beyond  which  they  do  not  increase,  fig2. 


-  2  - 


P]. 


In  the  present  work  the  nature  of  this  phenomenon  is  explained  in  terms  of  the  residual  thermal 
stresses  produced  by  cooldown  from  the  high  formation  temperature  of  over  1000  CO  to  room 
temperature  and  the  damage  in  form  of  microcracking  which  develops  during  load  cycling.  It  is 
shown  analytically  that  the  macro- residual  strain  can  be  expressed  in  terms  of  effective  thermal 
expansion  coefficients  (TEC)  of  the  composite  as  functions  of  the  damage  produced  at  the  vari¬ 
ous  stages  of  loading  history.  These  TEC  are  defined  in  the  secant  sense  for  the  temperature 
range  under  consideration.  The  theory  is  illustrated  by  evaluation  of  limit  residual  strains  for  uni¬ 
directional  fiber  composites  and  cross- ply  laminates. 


THEORETICAL  DEVELOPMENT 

When  a  homogeneous  stress-free  body  is  in  a  state  of  uniform  temperature  <f>o  and  the  tempera¬ 
ture  is  changed  to  another  uniform  state  the  body  remains  in  a  stress- free  state.  If,  however,  an 
initially  stress- free  composite  body  is  subjected  to  such  a  temperature  change,  internal  stresses 
develop  due  to  the  mismatch  of  elastic  and  thermal  expansion  properties  of  the  phases.  Such 
stresses  are  known  as  thermal  or  residual  stresses.  In  the  case  of  the  el^tic-brittle  composites 
considered  here  sufficiently  large  thermal  stresses  will  develop  damage  in  the  form  of  cracks, 
within  the  phase  regions  and/  or  on  the  interface.  Such  damage  produced  by  a  temperature  change 
alone  will  be  symbolically  denoted  as  ,  the  strain  and  stress  fields  at  this  temperature  will  be 
written  and  o<^(x,D^), respectively,  and  it  is  assumed  that  at  temperature  ^  there  are 

no  strains  and  stresses.  . 

A  thermoelastic  stress- strain  relation  for  temperature  dependent  properties  within  the  tempera¬ 
ture  range  [</io  ,</>]  may  be  expressed  in  the  alternate  forms 

ff  =  C(<t>)  e+T(<l>)  (<}>-4>o) 

€  =  S(</))  <r  +  a(«/>) 

r  =-C  a 

where  C,S  and  a  are  temperature  dependent  stiffness,  compliance  and  TEC  tensors.  It  is  implied 
that  a  is  a  secant  TET  which  is  defined  as  the  free  thermal  strain  produced  by  the  temperature 
change  .divided  by  the  temperature  change. 


(a) 

(b)  (1) 
(c) 
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An  effective  secant  TEC  a((f>)  of  a  composite  material  is  defined  in  analogy  with 
secant  a(<^)  in  terms  of  the  average  strain  produced  by  uniform  temperature  change  of  a  compo¬ 
site  which  is  not  loaded.  Spatially  uniform  temperature  is  a  basic  situation  for  composite  materials 
for  it  is  the  steady  state  temperature  of  a  composite  which  is  placed  in  a  uniform  temperature 
environment.  Since,  by  definition,  there  is  no  surface  load  the  average  stress  vanishes  and  the 
average  strain  for  temperature  change  which  may  be  called  the  macro- free  thermal  strain, 

is  written 


(2) 

where  the  TEC  is  also  a  function  of  the  eventual  damage  produced  by  thermal  stress.  It  is  seen 
from  (2)  that  the  computation  of  effective  TEC  requires  the  evaluation  of  a  strain  field  within  a 
cracked  composite  material  from  which  the  average  strain  can  then  be  determined.  The  form  of 
(2)  is  based  on  the  tacit  assumption  that  it  is  only  the  final  crack  geometry  which  is  of  impor¬ 
tance  and  not  the  history  of  crack  formation.  While  this  is  obvious  for  temperature  independent 
properties  it  may  not  be  evident  for  the  case  of  temperature  dependence  when  it  might  be  thought 
that  the  temperatures  at  which  the  various  cracks  appear  are  also  of  importance.  It  is  shown  in  the 
appendix  that  indeed  the  temperature  history  of  crack  development  is  of  no  importance.  One 
might  visualize  the  situation  as  permission  to  translate  cracks  which  were  formed  at  various  tem¬ 
peratures  to  formation  at  one  same  temperature.  Thus  we  may  formulate  what  may  be  called  a 
temperature  translation  principle  for  crack  formation  :  the  strain/ stress  fields  within  a  cracked 
composite  material  with  temperature  dependent  properties,  produced  by  the  temperature  change 
4>-4>o  ,is  a  function  of  the  final  crack  configuration  at  temperature  and  not  of  the  temperature 
history  of  crack  formation. 

Suppose  now  that  the  composite  is  subjected  to  load  defined  by  the  homogeneous  trac¬ 
tion  boundary  condition 

T(S)  =  oOn(S)  (3) 

where  S  is  the  bounding  surface,  o®  is  a  constant  stress  tensor,  which  is  also  the  average  stress  in 
the  composite,  and  n  is  the  normal  to  S.  The  loading  (3)  produces  additional  damage  and 
additional  strain  and  stress  £T(x,D<^+DT)  and  <rT(x,D<^+DT),  respectively.  Then  average  strain 
is  added  to  (2)  in  consequence  of  the  loading  which  is  composed  of  two  parts,  one  reversible  and 
the  other  irreversible.  The  reversible  part  is  of  the  form 

jT  =s*(<#),D<^-1-DT)oO  (4) 

where  S*  is  the  effective  compliance  tensor  of  the  damaged  composite.  When  the  composite  is 
unloaded  until  tP  =  0  ,assuming  that  no  additional  damage  is  produced  during  the  unloading. 
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the  strain  (4)  will  reduce  to  zero.  The  irreversible  part  of  the  strain  added  by  the  load  is  due  to 
the  fact  that  additional  damage  has  been  produced  by  the  load.  Thus  after  unloading  the  com¬ 
posite  material  is  in  a  new  state  of  stress  and  strain  associated  with  the  damage  state  D<#’-1-DT.  It 
follows  from  the  temperature  translation  principle  formulated  above  that  this  new  state  of  strain 
is  the  same  as  in  a  specimen  in  which  this  damage  state  would  have  been  produced  after  a  temper¬ 
ature  change  without  the  loading/  unloading  cycle.  Therefore  in  analogy  with  (2)  the  aver¬ 

age  strain  after  unloading  can  be  expressed  in  the  form 

€  =  (5) 

Specimens  which  are  tested  at  room  temperature  have  been  previously  cooled  down  in  manufac¬ 
ture  and  therefore  the  strain  state  (2)  is  considered  as  a  zero  strain  state.  Consequently,  the  load- 
unload  cycle  produces  a  macro- residual  strain  which  is  the  difference  between  (5)  and  (2).  Thus 

cr  =  [a*(0,D^+DT)  -  a*(<l>,D<i>)](<h<h)  (6) 

Thus  the  average  residual  strain  is  determined  by  the  effective  TEC  at  the  damage  states  before 
and  after  the  load  cycle. 

Similarly,  if  N  load  cycles  of  any  amplitudes  are  applied  and  the  cumulative  damage 
due  to  these  cycles  is  denoted  then  on  the  basis  of  the  arguments  leading  to  the  derivation  of 
(6)  the  average  residual  strain  will  be 

€'■  =[«*(</>, D'^+DN)  -  a*(</),D<^)](</^(/>o)  O) 

It  should  be  noted  that  from  temperature  translation  the  actual  details  of  the  loading 
history  which  produces  a  certain  damage  state  are  of  no  importance.  If  the  same  damage  is  pro¬ 
duced  by  different  loading  histories  -  the  residual  strain  will  be  the  same,  whether  the  loading  is 
only  one  cycle  or  many  cycles  of  different  amplitudes. 

It  often  happens  that  a  brittle  composite  material  becomes  crack  saturated.  This  implies 
that  for  increasing  number  of  load  cycles  it  becomes  more  and  more  difficult  to  produce  addi¬ 
tional  cracks  and  the  effective  TEC  approaches  a  limit  asymptotically.  An  example  of  this  for 
cracked  laminates  has  been  given  in  Hashin  (1988)  where  it  has  been  shown  that  the  effective 
TEC  of  cross-ply  laminates  approaches  asymptotic  limits  with  increasing  crack  densities  in  the 
layers.  It  follows  from  (7)  that  in  this  situation  the  residual  average  strjun  approaches  a  limit  . 
This  phenomenon  is  well  known  and  examples  based  on  actual  experiments  are  shown  in  figs.  1,2 

Note  that  in  accordance  with  the  discussion  above  the  value  of  also  depends  solely 
on  the  end  state  of  crack  saturation  and  not  on  the  history  of  formation. 
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APPLICATIONS 

It  follows  from  the  theory  developed  above  that  analytical  evaluation  of  macro- resi¬ 
dual  strains  requires  :  (a)  Knowledge  of  the  damage  produced,  (b)  Analytical  evaluation  of  the 
effective  TEC  as  a  function  of  the  damage.  Damage  can  be  described  on  the  basis  of  suitable 
observation  of  a  specimen.  It  would  of  course  be  very  desirable  to  be  able  to  predict  damage  ana¬ 
lytically.  Work  of  this  nature  has  been  described  in  Nairn  and  Hu  (1994)  and  Hashin  (1996). 
Since  .however,  this  subject  is  still  in  its  infancy  it  will  not  be  considered  here.  The  examples  to 
be  given  here  will  be  concerned  with  the  relatively  simple  case  when  the  the  temperature  change 
does  not  produce  damage  and  the  subsequent  load  cycling  is  continued  until  the  specimen  is  crack 
saturated.  Thus  the  purpose  is  to  evaluate  the  limiting  final  value  of  macro- residual  strain. 

As  a  simple  first  example  consider  a  one-dimensional  model  of  a  unidirectional  fiber 
composite,  fig.  3,  where  >of.  Here  m,f  denote  matrix  and  fiber,  respectively  and  the  TEC  are 
in  fiber  direction.  Assume  that  cooldown  from  <}>o  to  <l>  produces  no  damage.  Then  the  effective 

TEC  in  fiber  direction  is  given  with  good  approximatiom  by 


4(</.,D<^>)  =4(<#.,0)  = 


+  ojfEfVf 
+Ef  Vf 


(8) 


where  Vn,  and  Vf  are  the  volume  fractions  and  the  Young's  moduli  are  un  fiber  direction.  More 
accurate  expressions  are  available  in  the  literature  ;  e.g.  Rosen  and  Hashin  (1970),  Hashin 
(1979). 


The  composite  is  now  loaded  by  tensile  average  stress  (P  in  fiber  direction.  This  pro¬ 
duces  tensile  stress  in  the  matrix  which  is  superposed  on  the  tensile  matrix  stress  produced  by 
cooldown  (since  by  hypothesis  £Sn  >or)-  Assume  that  the  specimen  is  loaded  to  the  point  of 
matrix  crack  saturation.  This  matrix  crack  accumulation  is  the  damage  D^.  The  matrix  then 
becomes  ineffective  and  therefore  the  effective  TEC  reduces  to  Thus 

a;(<#.,D<^-fDT)  =  ci(<^,,DT)  =  Of  (9) 


It  follows  from  (7-9)  that 
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(10) 


(<f^<l>o) 


Since  <!>  <<f>o  and  ^  <qiij  this  is  a  positive  residual  strain. 

The  second  example  is  concerned  with  a  symmetric  cross-ply  laminate,  fig.4,  where 
the  layers  are  made  of  unidirectional  fiber  composites.  The 

laminate  is  cooled  down  from  stress-free  temperature  ^  to  temperature  If  no  damage  is  pro¬ 
duced  by  the  cooldown,  evaluation  of  the  laminate  TEC  and  internal  stresses  is  elementary.  For 
the  case  of  equal  layer  thicknesses  the  TEC  has  the  simple  form 


«A  +  (ar-«A) 


^+i'A 

1  H“Ey^/  E^p 


(11) 


where  ^*^1  (fiber  direction)  and  transverse  TEC  of  the  unidirectional  layer  mater¬ 

ial,  Ey^  and  Ej  are  the  Young's  moduli  for  these  directions  and  is  the  axial  Poisson's  ratio. 


Denoting  the  in-plane  axial  (in  fiber  direction)  and  transverse  layer  stresses  and 
Op,  respectively,  we  have 


_  Ea  («r-«A)  . . 

“  1+2»/a+Ea/Et 

Op  =  -Oa 


(12) 


All  TEC  in  (11,12)  are  to  be  interpreted  in  the  secant  sense  with  respect  to  tempera¬ 
ture  change  from  stress- free  temperature  0o  to  E^al  temperature  <j). 

In  unidirectional  fiber  composites  with  polymeric  or  metallic  matrix  generally  op  >aA 
.  It  follows  from  (12)  that  when  laminates  made  of  such  materials  are  cooled  down  from  4>q, 
tranverse  tensile  stresses  will  develop.  Such  stresses  may  crack  all  layers  along  fiber  directions 
producing  so-called  intralaminar  cracks.  In  order  to  keep  the  anlysis  simple  it  will  here  be 
assumed  that  cooldown  does  not  produce  cracks.  Next  the  laminate  is  subjected  to  tensile  cyclic 
uniaxial  stress  in  a  direction  transverse  to  the  fibers  in  the  inner  layer,  say,  which  is  designated 
as  the  X  direction.  This  load  produces  transverse  tensile  stresses  in  the  inner  layer  which  super¬ 
pose  on  the  tensile  thermal  stresses  due  to  cooldown.  It  is  assumed  that  cracks  are  produced  by 
this  cycling,  fig.  4,  and  as  is  known  from  experiment,  after  a  sufficient  number  of  cycles  crack 
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density  in  the  inner  layer  becomes  so  large  that  no  more  cracks  can  be  added  by  the  load  cycling 
which  implies  that  a  stage  of  crack  saturation  has  been  reached.  The  effective  TEC  in  load  direc¬ 
tion  for  this  limiting  case  has  been  given  in  Hashin  (1988)  in  the  form 


a*  =  (1  +  ^o) 

where 

Bq  =  -(1  +  1/X)  j'a/Ba 
Co  =  1/Ea  +  1/XBr 

(14) 

X  =  12/  tj 

kj  =<4x/<^ 


Here  inner  and  outer  layers  have  been  labeled  1  and  2,  respectively,  tj  and  t2  are  layer 
thicknesses  and  <4x  the  stress  in  the  inner  layer  for  an  undamaged  laminate.  Expression  (13) 
can  be  very  much  simplified  to  read 


ot  =  OfA 


Pa.  (ttr  -  «a) 
17 +X 


V  “  Ea  /  E j 


(15) 


Since  ot^  is  an  unwieldy  expression  for  arbitrary  X,  we  shall  restrict  ourselves  to  the  simple  case 
of  X=l.  Then,  identifying  (11)  with  a  (<^,D*)  and  (14)  with  X=1  as  a*(</.,D<^+DT)  ^^here  in 
this  case  there  is  no  thermal  damage,  we  have  from  (6)  after  some  algebra 


6** 

^XX 


(ot-«a)  (i+V-2pl)  , .  . 
--  {l+i])il+ri+2v^) 


(16) 


and  it  is  seen  that  (16)  is  a  positive  strain  increment  ,since  off  >aA  and  <l>  <<^>o. 

There  is  of  course  no  difficulty  to  carry  out  all  of  the  preceding  calculations  for  the 
case  of  unequal  thicknesses.  This  is  best  done  numerically  as  the  analytical  results  become 
unwieldy. 
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The  situation  op  <oc/^  may  arise  in  ceramic  composites  ;  for  example  :  the  TEC  of  a 
material  composed  of  Nicalon  fibers  and  SiC  matrix.  In  that  case  the  stresses  (12)  reverse  sign 
and  Op  is  a  compressive  residual  stress  which  can  not  produce  cracking.  If  the  laminate  is  now 
subjected  to  tensile  cycling  the  compressive  residual  transverse  stress  subtracts  from  the  tensile 
stress  produced  by  external  load.  The  applied  tensile  stress  may  be  large  enough  to  produce 
cracking  but  after  unloading  the  inner  layer  will  again  be  in  compression  and  the  cracks  will  be 
closed.  Therefore  there  is  no  additional  crack  effect  with  respect  to  the  state  before  loading  and 

there  will  thus  be  no  residual  macro- strain. 

Results  of  sample  calculations  of  limit  residual  strain  and  layer  properties  for  a  cer¬ 
amic  matrix  composite  are  shown  in  the  table  below.  The  material  consists  of  T300  graphite 
fibers  embedded  in  SiC  matrix.  The  volume  fraction  of  fibers  is  0.4  ,the  reference  temperature  is 
<t)Q  =  1025  CO  and  the  laminate  is  cooled  down  to  a  temperature  of  <^»  =  25  C®. 


Table 


Material 

Ea 

Et 

va 

“A 

«r 

OA 

oy 

T300/SiC 

434.5 

166.1 

0.179 

4.45 

7.34 

-25.54 

25.54 

7.14 

GPa 

GPa 

10-6/ CO 

10-6/ CO 

MPa 

MPa 

o 

t 

CONCLUSION 

It  has  been  shown  in  general  fashion  that  macro- residual  strains  produced  in  elas¬ 
tic-brittle  composites  under  cyclic  loading  are  tied  to  the  internal  thermal  residual  stresses  and  the 
damage  produced  by  the  cyclic  load  and  that,  furthermore,  these  residual  strains  can  be  analyti¬ 
cally  expressed  in  terms  of  effective  TEC  as  functions  of  damage.  The  problem  of  evaluation  of 
residual  strains  has  thus  been  reduced  to  the  well  known  problem  of  evaluation  of  effective  TEC. 

Some  simple  illustrative  evaluation  examples  have  been  given  and  there  is  no  diffi¬ 
culty  to  carry  out  more  complicated  calculations  for  cross- ply  laminates  in  which  all  plies  are 
cracked.  The  theory  can  also  be  readily  applied  numerically  on  the  basis  of  numerical  calculation 
of  effective  TEC  for  cases  of  interest. 
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APPENDIX 


The  governing  differential  equations  of  any  homogeneous  constituent  phase  of  the 
temperature  dependent  thermoelastic  composite  body  are  obtained  by  introduction  of  (la)  into  the 
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stress  equilibrium  equations.  Thus 

CC^ijuUM  =  0  (Al) 

where  u(x)  is  the  displacement  and  there  is  no  contribution  from  the  thermal  term  in  (la)  since  it 

is  constant  and  thus  its  derivatives  vanish. 

From  (la)  the  traction  vector  associated  with  surface  normal  n  may  be  written 

Ti  =  CjjkiUk.inj  +  r(<^) 

The  traction  (A2)  and  the  displacement  u  must  be  continuous  at  the  phase  interfaces. 

Let  it  be  assumed  that  the  composite  is  subjected  to  any  temperature  variation 
within  the  range  [<^,<#>]  and  a  surface  traction  distribution  T(S,<^).  It  follows  from  the  mathemati¬ 
cal  nature  of  the  composite  boundary  value  problem  (BVP)  described  above  that :  (a)  The  solu¬ 
tion  of  the  BVP  for  any  temperature  is  a  function  of  ^  and  not  of  temperature  history,  (b)  Since 
the  BVP  is  linear  in  space  its  solution  u,<rat  any  temperature  <l>  can  be  constructed  as  the  superpo¬ 
sition 


u  =  u^  +  u*^ 


(A3) 


O'  = 

where  the  <l>  superscript  indicates  the  solution  of  the  BVP  for  the  temperature  input  with  T(S)  = 
0  and  the  T  superscript  indicates  solution  for  specified  <t>  with  load  input 

We  now  describe  two  different  scenarios  for  the  same  composite  body  with 
bounding  surface  S.  First,  the  temperature  is  changed  from  to  <^>2  ,  and  a  crack  with  surfa^  Sc 
appears  at  temperature  Second,  the  temperature  is  changed  from  to  the  same  identical 
crack  at  same  location  appears  at  temperature  <^i  and  then  the  temperature  is  changed  to  We 
wish  to  show  that  in  the  two  cases  the  internal  fields  at  temperature  ^  are  identical. 

Consider  first  the  case  when  the  composite  body  undergoes  a  temperature  change 
ij)Q  -*  (j)2  ,  there  is  no  load  on  S  ,  i.e.  T(S)  =  0,  and  no  crack  appears.  The  stress  field  at  any  tem¬ 
perature  <l>  for  this  case  is  identified  with  o^(<^)  in  (A3).  The  traction  on  the  anticipated  crack 
surface  Sc  which  at  this  stage  is  a  mathematical  surface  is  denoted  T(Sc,0).  Next  we  consider  the 
same  body  with  no  load  on  S  and  no  temperature  change  which  is  subjected  to  the  traction 
-T(Sc,<^>)  on  Sc.  The  stress  field  produced  by  this  load  is  identified  with  <rT(<^)  in  (A3).  The 
stress  field  of  the  first  scenario  is 
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The  stress  field  of  the  second  scenario  is 

+  <^(^)~  <^(<^*1)  (^5) 

where  the  last  two  terms  in  (A5)  are  due  to  the  temperature  change  <j)i  -*  <t>2  .  In  each  case  the 
tractions  associated  with  the  stress  field  vanish  on  the  crack  surface  Sc  and  on  the  bounding  sur¬ 
face  S.  It  is  seen  that  the  two  scenarios  result  in  the  same  stress  field  and  similarly  in  the  same 
displacement  and  strain  field. 

Now  the  same  reasoning  can  be  applied  again  to  two  scenarios  of  an  additional 
crack  appearing  at  different  temperatures  to  conclude  finally  by  induction  that  temperature  trans¬ 
lation  of  cracks  is  valid  for  any  number  of  cracks. 


Stress  (MPa) 


TENSILE  STRESS  (MPa) 


TENSILE  STRAIN  (%) 


Fig.2  -  Limit  residual  strain  in  cyclic  loading  of  ceramic  matrix  composite 


Fig.3  -  Unidirectional  fiber  composite  with  cracked  matrix 
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APPLICATION  TO  LAMINATE  CRACKING  ANALYSIS 


Z.  Hashin 


Finite  Thermoelastic  Fracture  Criterion  with 
Application  to  Laminate  Cracking  Analysis 


Introduction 

The  dominant  mode  of  damage  development  in  composite  materials  with  polymeric  or 
ceramic  matrix  is  the  accumulation  of  cracks.  In  particulate  composites  and  in  unidirectional 
fiber  composites  such  cracks  occur  as  matrix  cracks,  interface  disbonds  and  fiber  breaks.  In 
fiber  compKJsite  laminates  so-called  intralaminar  cracks  accumulate  within  the  layers,  along 
the  fibers.  It  is  of  primary  importance  to  develop  an  anal3rtical  capability  for  prediction  of 
sudi  crack  accumulations  and  it  is  with  this  goal  that  the  present  work  is  concerned. 

The  subject  is  within  the  realm  of  fracture  mechanics  but  the  approach  is  of  necessity 
different  from  the  one  adopted  in  classical  fracture  mechanics.  The  latter  is  primarily  con¬ 
cerned  with  differential  extension  of  a  sin^e  crack,  or  perhaps  of  several  cracks  with  well 
defined  locations.  In  composites  the  concern  is  with  spontaneoxis  appearance  of  many  cracks 
of  finite  surface  and  ill  defined  locations.  It  may  therefore  be  appropriate  to  refer  to  the  sub¬ 
ject  as  finite  fracture  mechanics.  Spontaneous  crack  formation  implies  that  new  cracks 
appear  in  very  short  time  and  it  is  not  possible  or  of  interest  to  follow  the  history  of  their 
development.  The  formation  of  new  cracks  can  thus  be  regarded  as  fracture  events.  No 
single  crack  produces  failmre  and  every  cradc  is  arrested  by  a  boundary.  The  imcertainiy  of 
location  of  cracks  introduces  as  important  geometrical  parameter  the  crack  density  and 
perhaps,  as  will  be  seen  later,  other  global  quantities  such  as  standard  deviation  of  crack 
spacings. 

Furthermore,  composites  are  manufactured  at  relatively  high  temperatmre  and  are  sub¬ 
sequently  cooled  down  to  room  temperature.  This  process  introduces  residual  stresses, 
possibly  of  considerable  magnitude,  which  cannot  be  ignored  and  mxist  be  incorporated  into 
the  formulation  and  analysis. 

There  are  two  major  tasks  in  the  analysis  of  the  phenomena  described  above;  (a)  Es¬ 
tablishment  of  a  criterion  for  new  crack  formation,  (b)  Stress  analysis  of  the  cracked  body 
for  purpose  of  application  of  the  fracture  criterion.  A  cxirrent  major  research  activity  within 
the  present  subject  involving  these  two  tasks,  is  analysis  of  intralaminar  cracking  of  fiber 
composite  laminates  which  has  grown  into  a  voluminous  literature.  A  recent  review  has  been 
given  by  Naim  and  Shu  (1994). 

Literally  all  of  the  work  on  cracked  laminates  has  been  concerned  with  so-called  cross 
ply  laminates  which  consist  of  layers  with  reinforcement  directions  of  0®  or  90®.  The  most 
commonly  used  method  of  analysis  has  been  shear  lag  ,  which  appears  to  have  been  initi¬ 
ated  for  cracked  laminates  by  Garrett  and  Bailey  (1977).  However,  shear  lag  is  a  roughly 
approximate  one  dimensional  method  of  analysis  which  brings  in  an  artificial  shear  transfer 
zone  in-between  the  laminate  layers.  The  thickness  of  this  zone  is  an  vuiknown  parameter 
which  can  only  be  obtained  from  some  experimental  information.  Furthermore,  the  shear 
stresses  as  calculated  by  shear  lag  do  not  vanish  on  crack  surfaces,  as  they  should. 

A  two  dimensional  variational  method  of  stress  analysis  has  been  established  by  Hashin 
(1985)  based  on  the  sole  approximation  that  in-plane  layer  stresses  are  constant  over  layer 
thickness.  This  has  been  applied  for  laminates  with  one  kind  of  layers  cracked  and  has  been 
generalized  in  Hashin  (1987)  to  the  case  when  all  layers  are  cracked,  which  required  a  3- 
D  variational  analysis.  McCartney  (1992)  has  constmcted  approximate  elasticity  solutions 
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which  yield  resiilts  similar  to  the  ones  obtained  in  Hashin  (1985,1987).  Varna  and  Berglvmd 
(1994)  have  generalized  the  variational  method  by  permission  of  certain  variations  of  in¬ 
plane  stresses  over  layer  thickness,  at  the  price  of  considerable  complexity;  Naim  (1989)  has 
generalized  the  Hashin  (1985)  analysis  to  thermoelasticity. 

Regarding  fracture  analysis,  there  is  here  adopted  the  traditional  point  of  view  that  a 
fracture  criterion  should  be  based  on  energy  release.  According  to  Naim  and  Shu  (1994)  this 
approach  was  first  applied  to  laminate  cracking  by  Parvizi  et  al  (1978)  and  has  been  subse¬ 
quently  used  in  many  investigations  which  tiniversally  employed  the  shear  lag  approximation. 
Naim  (1989)  was  the  first  to  use  the  much  more  accurate  variational  method  in  conjunction 
with  energy  release  criteria  to  predict  crack  densities  in  cracked  laminates,  obtaining  good 
agreement  with  experimental  data;  Naim(1989),  Liu  and  Naim  (1992). 

In  the  present  work  it  is  first  shown  that  the  energy  release  required  for  formation  of 
new  finite  crack  surface,  by  load  and/or  temperature  input,  in  a  body  which  contains  an 
initial  field  of  residual  stress  can  be  expressed  entirely  by  stress  (or  strain),  thus  providing 
a  simple  and  convenient  fracture  criterion.  Furthermore,  on  the  basis  of  this  result,  it  is 
shown  that  an  upper  bound  on  the  energy  release  is  obtained  with  admissible  instead  of 
actual  stresses.  Then  a  new  simple  method  for  thermoelastic  analysis  of  cracked  laminates 
is  established  which  yields  thermoelastic  results  by  a  simple  replacement  in  corresponding 
isothermal  results.  On  the  basis  of  the  above  results  crack  formation  criteria  for  cross-ply 
laminates  are  established  in  terms  of  statistical  parameters  of  crack  distribution.  It  is  shown 
that  these  criteria  assume  simple  forms  for  small  and  large  cradc  densities. 


Thermoelastic  Cracking  Criterion 


This  section  will  be  concerned  with  the  development  of  a  criterion  for  the  spontaneous 
formation  of  new  cracks  due  to  smface  load  on  a  thermoelastic-brittle  body  in  which  there 
are  residual  stresses  due  to  previous  temperatrue  history.  Unlike  classical  fracture  mechanics 
where  infinitesimal  increase  of  crack  surface  is  considered  the  concern  is  here  with  finite 
increase. 

Consider  an  elastic-brittle  composite  material  body  which  is  first  subjected  to  temper¬ 
ature  change  T  without  load  and  is  subsequently  loaded  by  surface  tractions  T(S')  on  the 
hounding  surface  S.  For  composites  the  most  common  and  important  case  is  spatially  con¬ 
stant  temperature.  Because  of  thermal  expansion  mismatdi  the  temperature  change  will 
produce  residual  thermal  stress  while  the  subsequent  mechanical  load  will  add  mechan¬ 
ical  stress  resulting  in  an  initial  stress  field  a®  =  Here  and  from  now  on 

boldface  symbols  denote  tensors  and  vectors.  Initial  formation  of  internal  cracks  may  have 
occurred  due  to  temperature  change  and/or  load  and  their  surface  is  denoted  Sco-  Now 
consider  the  possibility  of  spontaneous  formation  of  additional  cradcs  for  same  surface 
load  T(iS}.  Then  the  thermal  and  mechanical  stress  fields  change  to  a’’  and  ,  the  final 
total  stress  field  is  =  O’*"  -f-  tr”*  and  the  final  crack  surface  is  denoted  Sc~  Denoting  by  T 
the  tractions  associated  with  a  and  all  differences  between  final  and  initial  quantities  by  the 
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(1) 


symbol  A,  (e.g.  Aa  =  a  -  a®;  ASc  =  Sc-  Sco),  there  are  the  following  relations 

T«(5co)  =  0;  T(5c)  =  0;  AT(5)=0 


The  stress-strain  relations  are 


€»•  =  So-’’  +  aT  ;  6’"  =  So-*" 

A€’’  =  SAo’’  ;  Ae’"  =  SA<r’' 


(2) 


where  S  is  the  elastic  compliance  tensor  and  a  is  the  thermal  expansion  tensor.  Frequent 
iise  will  be  made  of  the  theorem  of  virtual  work  for  a  cracked  body.  This  is  written  in  the 

form  f  r  f 

f  aedV  =  /  TudS'  -I-  /  TAidS  (3) 

Jv  Js  JSc 

where  Sc  is  the  cracks’  surface  considered  twice  as  if  the  cracks  were  very  fiat  voids,  a  and 
€  are  related  or  unrelated  equilibrated  stress  and  compatible  strain,  respectively,  and  u  is 
displacement.  The  contribution  from  the  crack  surface  integral  vanishes  in  the  following 
instances: 

W 

T  and  u  are  both  continuous  across  Sc 


Development  of  Fracture  Criterion 


"When  new  crack  surface  ASc  of  area  AAc  forms  without  change  of  surface  load  the 
internal  stress  and  strain  fields  change.  There  is  a  resulting  change  in  internal  mechanical 
energy  which  will  here  for  reasons  of  subsequent  development  be  expressed  in  terms  of  stress 
gngrgy  and  will  be  denoted  AU^-  Since,  by  hypothesis,  the  surface  load  does  not  change 
during  the  process  of  new  crack  formation  the  external  energy  change  is 

AU^  =  j^TAudS  (5) 

The  internal  mechanical  energy  change  will  be  defined  in  terms  of  stress  energy  with  density 
W,  defined  as 

W{a- )  =  ^  trSa  (6) 

Then  the  internal  energy  change  is 

ATT  =  AWdV  =  jT  [W{(t)  -  W(<T®)]  dV  (7) 

If  in  addition  energy  of  amount  AP  is  required  to  open  up  the  new  crack  surface,  then  the 
criterion  for  new  crack  formation  is  ^ven  by  the  energy  balance 

AU^  =  ALT  +  Ar  (8) 
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where,  1^  is  the  critical  energy  release.  It  may  be  assumed  that  there  exists  a  material 
surface  energy  per  imit  area  of  crack  opening  7,  in  which  case 

Ar  =  7AA‘^  (9) 


It  follows  that  7  is  also  the  analogue  of  the  critical  energy  release  rate  which  is  usually 
defined  as  the  derivative  of  the  energy  release  with  respect  to  crack  area. 

From  (1,  3,  4,  5) 

AJ7^=  /  o-Aedy  (10) 

Jv 

and  using  (1)  this  can  be  expressed  as 

=  jT  (o-»SAo-  +  AaSAa)dV  (11) 

The  internal  energy  (7)  can  be  expanded  to  read 

AU^  =  jT  (<r“SA<T  +  ^A(TSAtr)dV  (12) 

It  then  follows  from  (8,  11,  12)  that 

Ar  =  -  f  AaSAadV  (13) 

2  Jv 

Thus,  the  critical  energy  release  has  been  expressed  entirely  in  terms  of  stress. 

In  view  of  (7, 12),  e<pi.  (13)  can  be  written 

Ar  =  jT  (aw  -  <7°SA<r)  dV  (14) 

The  second  integral  in  (13)  can  be  transformed  as  follows 

J  (r°SA<rdV  =  AaSa^dV  =  (e®  -  otT)  dV 


Now  from  (1  -  4) 

AaeW  =  jT  ATviW  =  (T  -  T®)u«d5'  =  -  jT  T Vd5"  =  0 

Therefore  .  . 

/  a°SAcrdV  =  -  Aa  otTdV 
Jv  Jv 

and  an  alternative  form  of  AIT^  is 

Ar  =  f  [W{a)  -  W(<rO)  +  (<t  -  o-°)aT]  dV 

J  V 

There  is  an  important  consequence  of  the  representation  (16).  If  the  surface  load  T(5')  is 
achieved  in  a  monotonic  sequence  of  loads  T\  T^,  •  •  •  ,  T",  •  •  •  ,  T,  at  the  same  temperature, 
and  if  (9)  is  valid  with  fixed  7,  then  for  any  load  increase,  T”"^  to  T”,  say 


(15) 


(16) 


AI^  =  7A4,  =  -  W(<t"-^)  +  (o-”  -  <7"-i)ar]  dV 


(17) 
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Siimmation  of  (17)  over  all  sequential  loading  steps  recovers  the  result(16)  but  with  an 
important  difference  of  interpretation.  Previously  the  stress  fields  <7°  and  a  where  associated 
with  same  surface  tractions  while  in  the  present  case  <7®  is  associated  with  surface  load 
T*(5)  ,say,  which  for  example  may  initiate  cracking,  while  <7  is  associated  with  T{S)  and 
a  subsequent  crack  geometry  produced  during  the  sequential  cracking  process.. 

There  is  of  course  no  difficulty  to  express  (13,  16)  in  terms  of  strain  by  use  of  the  relations 

(2). 

In  summary,  if  surface  load  T(5)  admits  increase  of  crack  surface  ASc,  then  the  stress 
increments  due  to  increase  of  crack  surface  must  satisfy  (13,  16)  where  is  the  critical 
energy  release  required  to  produce  the  new  crack  surface.  The  criterion  includes  as  special 
cases  the  situation  when  there  are  no  residual  stresses  and  on  the  other  hand,  the  situation 
when  cracking  is  due  entirely  to  residual  thermal  stresses  and  there  is  no  mediamcal  load.  In 
the  latter  case  (13,  16)  becomes  a  criterion  for  crack  formation  due  to  temperature  change 
alone.  When  energy  required  to  open  cracks  is  proportional  to  cradc  surface  area,  as  ex¬ 
pressed  by  (9),  the  criterion  (13,16)  is  also  valid  for  the  case  when  the  total  new  crack  area 
opens  up  sequentially  during  the  loading  process. 

It  should  be  further  noted  that:  (a)  Kinetic  energy  associated  with  the  dynamics  of  spon¬ 
taneous  cradc  growth  has  been  ne^ected.  (b)  The  results  obtained  are  valid  for  temperature 
dependent  properties.  Thus,  if  T  is  the  temperature  change  relative  to  initial  temperature 
To,  then  the  compliance  S  should  be  taken  at  temperature  T -|-To  while  the  thermal  expansion 
tensor  oc  should  be  defined  in  the  secant  sense  over  the  interval  (To  To  -f  T).  (c)  Everything 
done  here  is  valid  for  residual  stresses  which  are  not  necessarily  the  result  of  temperature 
change,  but  it  is  necessary  that  their  increments  due  to  additional  crack  formation  obey  the 
relation  Ae''=  SAo’’’. 


Bounding  of  Energy  Release 


The  complementary  energy  of  an  elastic  body  with  internal  stress  field  <7®(x),  temperature 
change  T(x)  and  tractions  T  prescribed  on  its  entire  external  and  internal  surfaces  is  defined 

^  /•  1  1  T^ 

Vc{<t’‘)  =  +  aa>‘T  -  jCp— MK  (18) 

where  To  is  a  reference  temperature  and  Cp  is  the  specific  heat  at  constant  stress.  Let 
it  be  assumed  that  due  to  internal  changes,  e.g.  additional  crack  formation,  the  internal 
stress  field  changes  to  cr(x)  while  the  temperature  field  remains  unchanged.  Then  the  new 
complementary  energy  is  Uc  (o')-  Now  consider  instead  of  the  actual  a  an  admissible  stress 
field  <7®  which  satisfies  equilibrium,  traction  continuity  and  boundary  conditions,  but  not 
compatibility,  and  define 


=  Vc{(T^)  =  Uc{<T^)  -  Uc{tT  ») 

Ao-“  =  <7“  -  <7°  AT®  =  T“  -  T® 
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It  follows  from  (18,  19)  that) 

AU^  =  jT  A<r“  (S<t"  +  otT)dV+^J^  Aa^SAtrW  (20) 

The  parenthesis  in  the  first  integral  is  recognized  as  €°  and  therefore  this  integral  may  be 
transformed  by  virtual  work  into 

f  AtT^e^dV=  f  AT“u®dS'+  f  AT^u^dS  (21) 

Jv  Js  J  Sc 

The  integral  (21)  vanishes,  the  proof  being  entirely  analogous  to  the  one  given  for  the  van¬ 
ishing  of  the  second  integral  in  (11).  Therefore 

i  /  Aa“SAo- W  =  f  W{Aa‘^)dV  (22) 

2Jv  Jv 

This  is  a  generalization  of  a  similar  result  ©ven  in  Hashin  (1985)  for  the  isothermal  case. 


Now  define 


AUc  =  Uc{tT)  -  Uc{<T^) 
where  <t  is  the  actual  stress.  Then  as  a  special  case  of  (22) 

(23) 

AUc  =  -  f  AaSAadV 

2  Jv 

(24) 

By  the  thermoelastic  principle  of  minimmn  complementary  energy 

U^  =  U^  +  AU^>Uc  =  U^  +  AUc 


and  therefore 

AU^  >  AUc  (25) 

Recalling  now  the  fracture  process  under  fixed  surface  tractions  described  above,  the 
stress  field  for  fracture  surface  Sco  is  identified  as  cr®  and  the  stress  field  for  Sc  =  Sco  + 
as  a.  Then  it  follows  from  (13,24)  that 

AUc  =  at  ,  AU$  >  AT^  (26) 

Thus  (22)  is  an  upper  bound  on  the  energy  release  AIT^. 

It  should  be  noted  that  the  proof  given  above  is  for  the  case  when  all  cracks  appear  at 
same  surface  load.  The  proof  is  easily  extended  to  the  case  when  cracks  form  progressively, 
at  different  loads. 
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Thermoelastic  Analysis  of  Cracked  Laminates 


Application  of  the  theory  developed  above  to  cracking  of  fiber  composite  laminates  re¬ 
quires  stress  analysis  of  a  cracked  laminate.  Xhe  laminates  to  be  considered  here  consist  of 
layers  or  plies  which  are  unidirectional  fiber  composites.  Xhe  present  work  is  concerned  with 
the  simple  case  of  [0®/90p]a  laminates,  also  known  as  cross-ply  laminates.  Xhis  implies  that 
fiber  directions  are  orthogonal,  n  and  p  are  the  number  of  plies  in  each  layer  and  s  stands 
for  symmetric  which  means  that  the  midplane  is  a  plane  of  symmetry.  Xhe  manufacturing 
temperature  of  the  laminate  may  be  identified  as  Tq  and  the  room  temperature  is  Tq  -f-  T. 
Xhus  T  is  a  negative  cooldown  temperature  change.  Xhe  case  to  be  considered  is  initial 
cooldown  of  the  laminate  and  subsequent  loading  by  constant  in-plane  force  Nj,,  fig.l. 


Cooldown  of  such  a  laminate  produces  tensile  stresses  in  transverse  direction  to  the 
fibers  which  may  produce  cracks  along  the  fibers,  in  the  0®  and  90®  layers.  Such  crack 
distributions  have  been  considered  in  Hashin  (1987)  but  for  reasons  of  simplicity  will  not  be 
considered  here  as  it  would  complicate  the  analysis  considerably.  It  will  thus  be  assumed 
that  cracks  accumulate  only  with  mechanical  loading.  Consequently  such  cracks,  whidi  are 
called  intralaminar  cracks,  appear  only  in  the  90®  layer,  for  the  transverse  stress  (Xyy  in  the 
0®  layer  for  Na  load  is  very  small. 

Xhe  usual  assumption  is  made  that  the  layers  can  be  considered  as  amsotropic  homo¬ 
geneous  with  effective  properties  of  the  unidirectional  fiber  composite.  Even  with  this  sim¬ 
plification  rigorous  analysis  of  the  stress  fields  in  such  a  cracked  laminate  does  not  appear 
possible.  Much  of  the  work  done  has  been  based  on  the  shear  lag  approximation.  More  accu¬ 
racy  has  been  achieved  on  the  basis  of  variational  analysis  in  terms  of  admissible  stress.  Xhe 
first  such  analysis  for  the  present  kind  of  cracked  laminate  and  for  isothermal  condition  has 
been  given  in  Hashin  (1985).  Xhis  analysis  has  been  extended  to  loading  after  temperature 
change  in  Naim  (1989).  In  the  present  work  the  thermoelastic  analysis  is  reconsidered  and 
is  presented  in  very  simple  form  which  is  particularly  convenient  for  evaluation  of  the  energy 
release. 

Xhe  average  stress  apphed  to  the  laminate  is 

(T  =  Nj2h  (27) 

When  there  are  no  cracks  the  stresses  in  each  layer  are  constant,  except  for  narrow  edge 
boundary  strips,  and  the  only  surviving  stresses  in  the  present  case  are  and  or®j,.  For 
temperature  diange  T  and  load  (27)  the  stresses  in  the  different  layers  may  be  written 

0^1;  =  <Ti  =  ki(T  =  02  =  -f  r^T  (28) 
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where  the  k  and  r  coefficients  depend  on  laminate  geometry  and  material  properties  and  the 
labels  1,2  from  now  on  indicate  the  90°  and  0°  layers,  respectively.  The  admissible  stress 
field  is  chosen  as  plane  in  the  xz  plane  with  cTxx  stresses  constant  over  layer  thickness.  Thus 


=  o-i[l  - 

Kl  =  M^)] 


(29) 


where  (pi  and  ^2  are  unknown  functions.  Insertion  of  (29)  into  the  section  force  equilibrium 
condition  in  x  direction  yields  the  relation  <p2  =  -{ti/t2)<pi.  Identifying  <pi  as  <p  and  using 
the  equilibrium  equations,  interlaminar  traction  continuity  and  external  boimdary  conditions 
[for  details  see  Hashin  (1985)]  it  follows  that  the  admissible  stress  within  a  typical  region 
between  two  cracks,  fig.2,  is  given  as 

=  -ai<P{x)  l^afx  =  {<rif>^)<Pi^) 

=  (Ti<p'{x)z  =  {<Ti/X)<p'{x){h  —  z)  (30) 


Ao-“]  =  ai<p"{x){hti  -  z^)l2  A(jf^  =  MX)(p"{x)(h  -  zfj’l 

Here  A  =  *2/^1,  primes  denote  x  differentiation  and  conforming  to  previous  notation  the 
stresses  (30)  must  be  added  to  the  stress  field  <r°  in  the  uncracked  laminate  to  obtain  the 
admissible  stress  a®.  Note  that  the  xz  and  zz  stresses  are  zero  in  the  imcracked  laminate. 

Since  normal  and  shear  stress  must  vanish  on  crack  surfaces  x  =  ±a  the  function  (p  must 
satisfy  the  boimdary  conditions 

^(±a)  =  1  ;  (p'i'^o)  =  0  (31) 


As  in  Hashin  (1985),  an  optimal  function  (p  is  constructed  by  utilization  of  the  principle 
of  minimum  complementary  energy.  It  follows  firom  (19)  that  instead  of  minimizing  it  is 
sufficient  to  minimize  (22).  This  requires  the  stress  energy  densities  of  the  two  layers  with 
respect  to  the  x,y,z  coordinate  system  in  terms  of  the  plane  stresses  (30).  Both  layers  are 
transversely  isotropic  with  axis  of  symmetry  y  for  layer  1  and  x  for  layer  2.  Therefore 


2W^(cr)  =  (t^/Et  —  Et  +  <^lzl  Er  + 

2W^  (cr)  =  ctIxIEa  —  "^(^xx^^zz^aIEa  +  olzlEr  +  ^"izI^A 


(32) 


where  the  elastic  ply  properties  in  (32)  are:  Ea,  ^a-  sodal  Young’s  modulus  and  associated 
Poisson’s  ratio;  Et  ,  I'r  -  transverse  Young’s  modulus  and  associated  Poisson’s  ratio; 
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Ga  »  Gt  -  axial  and  transverse  shear  moduli.  Then  for  the  region  —a 
Q<z<h, 


AUS:  =  r  r  W\Aa‘^^)dz  +  r  W^{A<T^^)dz 

J—a  Jo  Jti 


dx 


(33) 


After  the  z  integration  is  carried  out  minimization  of  the  functional  (33)  becomes  a 
classical  problem  of  the  calculus  of  variations  which  results  in  an  Euler-Lagrange  differential 
equation  for  the  function  (f).  Now  it  is  evident  from  previous  development  that  the  variational 
solution  for  the  thermoelastic  and  isothermal  elastic  cases  are  mathematically  identical. 
Indeed  the  only  difference  between  the  two  is  the  detailed  form  of  (28)  where  in  the  first  case 
r  ^  0  and  in  the  second  T  =  0.  Thus  the  isothermal  solution  developed  in  Hashin  (1985) 
can  immediately  be  used  to  write  down  the  thermoelastic  solution  simply  by  interpreting  ai 
as  (28).  But  before  this  is  done  it  is  necessary  to  take  into  account  that  intercrack  distances 
are  of  different  sizes.  Let  it  be  assumed  that  the  90®  layer,  of  length  L  m.  x  direction,  is 
cracked  into  N  pieces  where  any  intercrack  distance  is  denoted  2an.  Then  the  admissible 
Ao‘°  for  the  region  of  such  length  and  height  h  is  defined  by  (30)  with  boxmdary  conditions 


^„(±a)  =  1,  if„(±a)  =  0 


(34) 


According  to  Hashin  (1985)  the  Euler-Lagrange  equation  for  <^n  and  associated  boimdary 
conditions  are 


^4  d"  P  ^2  ^n{^pn)  —  1,  4^n{^Pn)  ~  0  (35) 

where 

^  =  x/ti,  pn  =  a^/h,  p  =  (Co2  -  C'ii)/C22,  q  =  G00/C22  (36) 

Coo  =  1/Et  +  1/XEa  C02  =  ^  (A  +  5)  - 

C^22  ==  (A  +  1)(3A^  +  12A  +  8)/60£t  Cii  = —(l/Gr  +  1/AG^)  (37) 

o 

This  determines  the  functions  <1)^  and  thus  the  admissible  stresses.  For  present  piurpose 
it  is  the  quantity  AUq  which  is  of  interest.  As  in  Hashin  (1985)  this  is  given  by 

ACJ  =  f;  AUi„  =  olilC22  Ex(p.)  (38) 

n=l  n— 1 

where 

xio  =  -no  (39) 

The  characteristic  equation  of  (35)  is 

r^+pr‘^  +  q  —  0  (40) 
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The  form  of  x  depends  on  the  nature  of  the  roots  of  (40).  When  these  roots  are  complex,  of 
the  form  ±(a:  +  ifi),  where  i  =  y/^  ,  then 


Xn  =  Xn  =  2Q:^(<a:^  +  0^) 


Cosh{2apn)  —  Cos{2/3pn) 
aSin{2Ppn)  +  j3Sinh{2apn) 


(41) 


If  the  roots  are  real,  thus  of  the  form  ±o:,  ±/?  ,then 


_  r  ^ _ 

Xn  -  Xn  -  Coth{apn)/a  -  Coth{0Pn)/0 


(42) 


Thus  the  form  of  the  function  Xn  depends  entirely  on  the  material  properties  and  stacking 
geometry  of  the  layers. 


Analysis  of  Laminate  Cracking 


When  the  laminate  described  above  develops  cracks  over  length  L  in  x  direction,  and  all 
cracks  appear  simultaneously  at  external  stress  a,  then  the  energy  release  per  unit  length  in  y 
direction  needed  to  open  up  these  cracks  is  bounded  from  above  by  (38).  In  view  of  previous 
ejqjerience  with  cracked  laminate  energy  calculation,  Hashin  (1985),  it  is  expected  that  (38) 
is  an  accurate  approximation  for  the  energy  release.  It  is  now  assumed  that  there  is  a  surface 
energy  7  per  unit  area  along  fibers,  which  is  a  material  property  of  the  imidirectional  ply 
material.  Since  by  hypothesis  the  layer  is  broken  into  N  pieces  there  are  N  cradcs  and 
therefore  Aj4c  =  {N  —  l)ti,  per  unit  length  in  y  direction  and  0  <  z  <  h.  Then  it  follows 
from  (9,  38) 

<^lhCii'Zx{l>^)  =  (N-l)l  (43) 

n=l 

If,  however,  cracks  initiate  at  some  external  average  stress  cTj  and  develop  progressively 
as  the  external  stress  increases  to  cr  then  it  may  be  shown  on  the  basis  of  the  interpretation 
of  17  for  progressive  cracking  that  43  is  modified  to 

[{a^  -  Oi)/E*  +  2{a  -  ai)alT  ]hJ2pn  +  olhCvi  X)  x{pn)  =  {N  -  1)7  (44) 

n=l  n— 1 

Since  N  is  generally  a  large  number,  N  —  1  may  be  replaced  by  N.  Then  (44)  can  be 
written 


[(<7^  -  afysz  +  2(a  -  c,)a:T  ](1  +  A)«,  P  +(^,,7  +  riT)\C^  X=  7  (45) 

where  an  overbar  denotes  mean  value  and  E*  ^  a*  are  the  effective  Yoimg’s  modulus  and 
thermal  expansion  coefficient,  respectively,  of  the  uncracked  laminate  in  x  direction. 
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The  numbers  p„  can  be  regarded  as  the  values  of  a  random  variable  p  with  probabihty 
density  function,  (PDF),  P{p).  Then  x  can  be  expressed  as 

X  =  P{p)x{p)dp  (46) 

The  specific  form  of  the  PDF  is  of  course  not  known  and  it  would  appear  that  the  choice 
should  be  based  on  experimental  data.  It  may  be  assumed  that  the  dioice  can  be  restricted 
to  those  PDF  which  are  specified  by  the  first  two  moments  of  p,  namely  the  mean  p  and  the 
standard  deviation  p,  and  in  that  case  X=X  (P,p)-  The  most  commonly  used  geometrical 
parameter  for  cradc  distribution  is  the  crack  density  c  which  is  the  number  of  cracks  per  unit 
length,  e.g.  mm.  Since  the  thickness  of  the  cracked  layer  is  2b\  it  follows  that 

p  =  \/2cti  (47) 

Thus  (45)  is  a  relation  between  the  input  stress  and  temperatiure  which  produce  the  cracking, 
the  geometrical  parameters  of  the  crack  distribution  and  the  surface  energy. 

It  is  instructive  and  helpful  to  consider  the  extreme  cases  of  small  and  large  crack  density 
which  are  characterized  by  large  and  small  p,  respectively.  A  typical  plot  of  the  function 
x{p)  for  a  graphite/epoxy  laminate,  where  the  relevant  form  of  x  is  (41),  is  shown  in  fig.  3. 
The  unidirectional  material  is  T300  graphite  fibers  in  fiberite  934  epoxy  with  fiber  volume 
fraction  0.55.  The  properties  of  the  material  are  :  =  128GPa,  Et  =  8.44GPa,  Ga  = 

3.85GPa,  i/x  =  0*3,  i^  =  0.35,  ax  = —0.55  10“®/^ >  =  35.9  10~®/C®.  It  is  seen  that 

the  plot  consists  of  a  linear  portion  for  a  range  0  <  p  <  pi,  a  curved  part  and  a  horizontal 
as3miptote  for  a  range  p2  <  p  <  oo.  Examination  of  such  plots  for  a  number  of  laminates 
made  of  fiber  composites  with  typical  fiber  volume  fractions  0.5  —  0.6  has  shown  that  in  aU 
cases  Pi  ~  0.9;  P2  ~  2.2.  This  also  seems  to  be  the  case  when  x  is  of  the  form  (42).  The 
nature  of  the  plot  is  illuminated  by  series  expansions  of  x  P  =  0  and  asymptotic  forms 
for  large  p.  Thus  for  small  p 

xS  =  (c?  +  0-‘f\p  -  +  +  0(/)J 

=  +  (48) 

And  for  large  p 

Xf  =  2a(a^  +  P^){1  -  2e"2"'’[Cos(2/3p)  +  ^Sin{2pp)\) 

xT  S  a/3(a  +  ^l)(l  -  (49) 
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The  absence  of  a  term  in  (48)  accoiints  for  the  substantial  initial  linear  part  of  the 
plot  and  the  rapid  decrease  of  the  escponentials  in  (49)  accounts  for  the  asymptote. 

The  asymptotic  part  of  the  x  curve  is  associated  with  crack  initiation  at  stress  a  = 

In  that  case  the  term  multiplying  P  in  (45)  vanishes  and  since  the  values  of  p  are  large  the 
negative  exponentials  in  (49)  are  neglected. 

This  leads  to  the  small  crack  density  or  crack  initiation  relations 


{hcTi  +  riTf2tiC22(xia^  +  i^)  =  7  X  =  X' 

(fci<Ti  +  riT)HiC^a^{a  +  /?)  =  7  X  =  X*” 


(50) 


According  to  previous  discussion  these  relations  should  be  valid  when  all  or  most  p„ 
are  larger  than  2.2.  The  actual  magnitude  of  small  crack  density  is  independent  of  crack 
producing  input  Oi.T,  Equ.  (50)  may  be  used  to  determine  in-situ  7  when  and  T  are 
known. 

It  should  be  noted  that  the  value  p  2.2,  when  x  becomes  flat,  fig.3,  does  not  imply 
small  cradc  density  for  the  in-plane  stress  When  the  cradcs  are  far  apart  this  stress 
builds  up  from  zero  on  the  crack  surface  to  the  value  <J^,  (28),  in  the  imcracked  laminate. 
Stress  variation  for  p  =  2  has  been  given  in  Hashin  (1985)  and  it  is  seen  that  in  this  case 
there  is  strong  crack  interaction  effect  and  builds  up  only  to  about  0.3<r^  midway  between 
cracks.  The  explanation  is  that  cr^  depends  on  <f>,  see  (29),  while  x  =  Therefore  when 

X  is  constant  (j)  is  cubic. 


Discussion  and  Conclusion 


Previous  approaches,  Naim  (1989)  in  terms  of  the  2-D  variational  solution.  Laws  and 
Dvorak  (1988)  in  terms  of  shear-lag  solution,  are  different  in  the  sense  that  it  has  been 
attempted  to  consider  the  details  of  progressive  development  of  the  crack  patterns  in  order 
to  evaluate  at  successive  stages  the  energy  releases  needed  to  produce  further  cracking.  This 
may  thus  be  termed  a  local  approach.  The  fundamental  question  which  arises  is:  given  the 
stress  variation  between  two  cracks  which  are  distance  2a„  apart.  What  is  the  probability 
of  appearance  of  a  new  crack  at  any  in-between  location  ?  But  in  reference  to  such  a 
question  it  is  necessary  to  consider  the  micro-mechanism  of  transverse  failure  of  a  ply.  There 
is  a  distribution  of  flaws  in  the  material  which  can  take  the  form  of  small  matrix  cracks, 
pores  and  fiber/matrix  interface  disbonds.  The  weakness  of  a  flaw  may  be  defined  by  the 
magnitude  of  the  stress  whidi  triggers  it  to  develop  an  intralaminar  crack.  The  weaker  the 
flaw  the  less  stress  is  required  to  initiate  a  nmning  crack  at  the  flaw  site.  Flaws  of  different 
weakness  are  randomly  distributed  along  the  layer  (  x  direction)  while,  by  comparison,  the 
stress  distribution  between  any  two  cracks  may  be  assumed  known  and  deterministic.  Even 
with  the  assumption  that  only  the  stress  and  not  the  other  components  are  involved 
in  the  failure,  the  problem  of  local  probability  of  failure  seems  intractable  as  nothing  is 
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known  about  the  statistics  of  flaw  distribution.  Wang  (1984)  has  considered  effect  of  flaw 
distributions  by  describing  them  as  microcracks  in  the  transverse  yz  plane,  of  size  much 
smaller  than  layer  thickness,  assuming  normal  PDF  for  crack  size  and  spacing.  But  those 
flaws  were  of  necessity  considered  as  classical  cracks  within  the  effective  fiber  composite,  in 
spite  of  the  fact  that  their  size  may  be  only  that  of  a  single  fiber,  which  negates  the  effective 
property  representation  of  their  surroimdings. 

The  most  simplistic  past  approach  has  been  semi-deterministic  in  that  it  has  been  as¬ 
sumed  that  a  new  crack  will  appear  midway  between  cracks  since  this  is  the  location  of 
largest  <J^.  This  approach  is  of  course  unacceptable  as  it  completely  ignores  the  presence 
of  flaws.  Naim  (1989)  and  Laws  and  Dvorak  (1988)  have  adopted  the  more  reasonable  as¬ 
sumption  that  the  probability  of  failure  at  any  location  x  is  proportional  to  the  magnitude 
of  However,  the  latter  analysis  is  based  on  the  unrealistic  stress  field  of  the  shear 

lag  approximation  and  therefore  the  results  obtained  by  Naim  should  be  much  preferred 
since  they  are  based  on  the  much  more  realistic  stress  field  as  determined  by  the  variational 
method. 

The  present  approach  is  related  to  that  of  Naim  (1989)  in  that  it  is  based  on  almost 
the  same  stress  analysis  of  a  cracked  laminate  but  there  are  important  differences.  The 
approach  may  be  termed  global  rather  than  local  since  no  attempt  has  been  made  to  follow 
the  complex  process  of  development  of  cradcs.  It  is  only  in  the  intermediate  range  of  curved 
X  plot  where  crack  geometry  information  beyond  crack  density  c  or  equivalently  p  is  needed, 
and  for  this  range  it  has  been  proposed  to  bring  in  the  required  additional  information 
through  the  standard  deviation  of  normalized  intercrack  distance  p. 

Another  distinguishing  featme  of  the  method  developed  here  is  the  expression  of  the 
fracture  criterion  (8)  in  terms  of  stress  only,  (13,  16).  This  establishes,  as  has  been  shown, 
that  the  expressions  (45,  50)  obtained  for  7  are  rigorous  upper  bormds  on  this  quantity.  In 
all  of  the  relevant  previous  literatiue  (5),  (8),  (9)  have  been  used  as  the  fracture  condition. 
This  requires  the  evaluation  of  displacements  from  an  approximate  stress  solution,  thus  from 
incompatible  strains,  and  thus  introduces  additional  approximation  of  uncertain  magnitude. 

The  usually  available  information  in  the  laminate  cracking  process  is  the  input  load  (x 
and  the  crack  density  c  which  determines  p.  There  is  of  course  no  fundamental  difficulty  to 
determine  p  from  section  photographs.  The  input  temperature  change  T  is  usually  taken 
as  the  difference  between  room  temperature  and  the  manufacturing  temperatiue  Tq  of  the 
laminate.  This  is  a  problematic  assumption  for  at  To  the  laminate  matrix  is  very  soft  and 
may  be  viscoelastic  rather  than  elastic.  In  view  of  the  elastic  behavior  requirement  of  the 
analysis  it  is  better  to  define  To  as  the  temperature  at  which  during  the  cooldown  process 
the  laminate  becomes  stiff  enou^  to  be  considered  elastic.  This  temperature  is  not  easy 
to  determine  directly  and  it  would  seem  that  it  should  be  backed  out  by  comparison  of  the 
analytical  results  derived  to  experimental  data  of  crack  density  versus  applied  load. 

The  most  important  physical  property  which  must  be  determined  is  the  surface  energy  per 
unit  area  or  critical  energy  release  rate  7.  The  simplest  possibility  is  to  regard  7  as  a  material 
constant  for  any  specific  unidirectional  fiber  composite  layer,  independent  of  layer  thickness. 
If  so  for  a  substantial  range  of  crack  density,  about  p  >  2.2  the  crack  density  is  independent 
of  input  a.  Thus  a  plot  of  crack  density  (ordinate)  versus  a  (abscissa)  shotJd  start  with  a 
vertical  straight  line  for  small  to  medium  crack  density.  Ebctensive  experiments  of  laminate 
crack  accmnulation  for  a  series  of  different  laminates,  made  of  different  materials  have  been 
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performed  by  Liu  and  Naim  (1992),  and  previously  by  Highsmith  and  Reifsnider  (1982) 
and  Wang  (1984).  There  are  cases  when  initial  crack  density  growth  is  load  independent  as 
described  above  but  in  other  cases  initial  crack  density  growth  is  accompanied  by  definite 
increase  in  a.  A  possible  explanation  for  the  latter  phenomenon  is  again  the  distribution  of 
flaws  in  the  material.  As  tr  is  applied  the  more  severe  flaws  wiU  be  activated  flrst  but  higher 
stress  will  be  required  to  activate  less  severe  flaws  to  result  in  cracks. 

The  theory  developed  should  be  compared  to  available  experimental  data.  For  this  pur¬ 
pose  it  is  necessary  to  have  crack  density  versus  applied  stress  data  for  several  different 
laminates  of  different  layup  geometries  (staddng  sequence)  made  of  the  same  layer  unidi¬ 
rectional  material.  The  values  of  7  and  T  must  be  backed  out  from  sudi  data  with  aid 
of  (50)  for  small  crack  density.  Then  verification  consists  of  ability  to  represent  the  other 
data  which  have  not  been  used  for  7,T  determination  with  these  values.  Such  verification 
has  been  given  by  Liu  and  Naim  (1992)  with  reasonably  good  results  on  the  basis  of  their 
analysis.  Fig.  4  shows  crack  densities  versus  load  for  3  different  laminates,  as  reported  by 
liu  and  Naim  (1992),  where  the  unidirectional  material  is  composed  of  T300  graphite  fibers 
and  Fiberite  934  epoxy  matrix.  The  values  of  p  ,  for  the  different  laminates  corresponding  to 
crack  density  c,  (47),  are  shown  near  the  experimental  data  of  each  laminate.  It  is  seen  that 
for  the  [0/90]a  and  [0/902]*  laminates  the  experimental  data  are  within  the  asymptotic  flat 
range  of  the  x  curve  when  crack  density  is  independent  of  the  value  of  applied  stress,  (50), 
and  the  experimental  data  verify  this  ,  with  some  scatter.  On  the  other  hand  the  [O/OOjs 
data  enter  into  the  range  of  curved  x  when  p  <  2.2  and  it  is  seen  that  these  data  deviate 
substantially  from  a  vertical  straight  line.  The  calculations  where  performed  with  the  values 
T  =  —160°,  7  =  1000  Joule/m.  Detailed  verification  of  experimental  results,  in  particular 
for  progressive  cracking,  is  deferred  to  the  futmre. 


Acknowledgments 


Support  of  AFOSR  through  the  Materials  Sciences  Corp.,  Fort  Washington,  PA,  Dr. 
Walter  Jones  contract  monitor,  and  the  Nathan  Cummings  Chair  of  Mechanics  of  Solids,  Tel 
Aviv  University,  is  gratefully  acknowledged.  Also  acknowledged  are  helpful  scientific  inter¬ 
changes  with  Professors  J.  A.  Naim  and  Y.  Benveniste  and  constmctive  critical  comments 
of  Professor  L.  Slepyan  and  Mr.  V.  Vinogradov. 


15 


REFERENCES 


Garrett,  K.W.  and  Bailey,  J.E.  (1977)  Multiple  transverse  fracture  in  90"  cross-ply 
laminates  of  glass-fibre  reinforced  polyester.  J.  Mater.  ScL  12,157-168. 

Hashin,  Z.  (1985)  Analysis  of  cradced  laminates:  A  variational  approach.  Mech. 
Mater.  4, 121-136. 

Hashin,  Z.  (1987)  Analysis  of  orthogonally  cracked  laminates  imder  tension.  J.  Appl. 
Mech.54,  872-879. 

Highsmith,  A.L.  and  Reifsnider,  K.L.  (1982)  Stiffness  reduction  mechanisms  in  com¬ 
posite  laminates.  Damage  in  Composite  Materials,  (ed.  Reifsnider,  K.L.)  pp. 
103-117.  ASTM  STP  115,  American  Society  for  Testing  Materials. 

Laws,  N.  and  Dvorak,  G.  J.  (1988)  Progressive  transverse  cracking  in  composite  lam¬ 
inates.  J.  Comp.  Mater.  22,  900-916. 

Liu,  S.  and  Naim,  J.A.  (1992)  The  formation  and  propagation  of  matrix  microcracks 
in  cross-ply  laminates  during  static  loachng.  J.  Reinforced  Plast.  Comp.  11, 
158-178. 

McCartney,  L.N.  (1992)  A  theory  of  stress  transfer  in  a  0"-90"-0"  cross-ply  laminate 
containing  a  parallel  array  of  transverse  cracks.  J.  Mech.  Phys.  Solids  40,  27-68. 

Naim,  J.A.  (1989)  The  strain  energy  release  rate  of  composite  microcraddng  :  A 
variational  approach.  J.  Comp.  Mater.  23,  1106-1129. 

Naim,  J.A.  and  Shu,  S.  (1994)  Matrix  microcracking.  Damage  Mechanics  of  Com¬ 
posite  Materials^  (ed  Talreja,  R.),  pp.  187-243.  Elsevier  Science  BN. 

Parvizi,  A.,  Garrett,  K.W.  and  Bailey,  J.E.  (1978)  Constrained  cracking  in  glass  fibre 
reinforced  eqoxy  cross-ply  laminates.  J.  Mater.  Sci.  13,  195-201. 

Wang,  A.S.D.  (1984)  Fracture  mechanics  of  sublamnate  cracks  in  composite  materials. 
Comp.  Technol.  Rev.  6,  45-62. 

Varna,  J.  and  Berglund,  L.A.  (1994)  Thermo-elastic  properties  of  composite  laminates 
with  transverse  cracks.  J.  Comp.  Tech.  Res.  16,  77-87 


16 


Fiberite  934/T300 


Fig.  4  -  Crack  Density  versus  Applied  Stress. 

Experimental  Data  from  Liu  and  Nairn 


INTERFACE  DEBOND  ANALYSIS 


Z.  Hashin  and  V.  Vinogradov 


Interface  Debond  Analysis 

Introduction 

Interface  debonding  in  composite  materials  is  a  problem  of  great  importance.  Formation 
of  interfacial  cracks  not  only  reduces  the  stiffness  of  composites,  but  also  produces  sig¬ 
nificant  local  stress  concentrations  which  may  result  in  appearance  of  microcracks  and, 
in  consequence,  produce  catastrophic  failure. 

In  this  connection  two  central  questions  are  under  discussion  in  the  present  chapter. 

1.  What  applied  mechanical  and  temperature  loads  produce  interface  debonding  in 
composite  materials  ? 

2.  What  parameters  of  a  composite  (such  as  volume  fraction  of  inclusions,  inclusion 
size  and  thermoelastic  properties  of  phases)  should  be  chosen  to  prevent  interface 
debonding  under  given  load  input  ? 

As  an  answer  to  these  questions  we  have  to  obtain  a  fracture  criterion  which  will  relate 
phases  and  interface  properties  and  other  details  of  the  composite  structure  to  critical 
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loads  and  temperature.  A  thermoelastc  fracture  criterion  developed  by  Hashin  [22]  is 
found  to  be  successfully  applicable  for  the  above-mentioned  purposes.  In  this  section  the 
method  will  be  briefly  outlined. 

Consider  a  brittle  elastic  composite  material  body,  subjected  to  an  external  mechani¬ 
cal  load  T(iS')  and  a  homogeneous  temperature  change  93,  which  may  be  applied  simulta¬ 
neously  or  sequentially.  Because  of  thermal  expansion  mismatch,  the  temperature  change 
produces  residual  stresses  in  the  composite  and  the  mechanical  load  leads  to  mechanical 
stresses  appearance  resulting  in  stresses  Initial  formation  of  internal  cracks  may  have 
occurred  due  to  temperature  change  and/or  loawi,  and  their  surface  is  denoted  Sco-  Now 
consider  the  possibility  of  spontaneous  formation  of  additional  cracks  for  the  same  surface 
load  T(5).  Then  the  stress  field  changes  to  a-  and  the  final  cr«w;k  surface  is  denoted  Sc- 
When  the  new  crack  surface  ASc  —  Sc  —  Sco  of  area  AAc  forms  without  change  of  the 
surface  load,  AT(5')  =  0,  the  external  energy  change  is 

AU^  =  jTAudS  (4.1) 

5 

The  internal  mechanical  energy  change  can  be  defined  in  terms  of  stress  energy  with 
density  W 

W{cr)——(rS<T  (4.2) 

The  internal  energy  change  can  be  written  as 

=  J  AW dV^  j [W(^)  -  W{a^)]dV.  (4.3) 

V  V 

If  an  additional  energy  of  amount  AF^  is  required  to  open  up  the  new  crack  surface,  then 
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the  criterion  for  new  crack  formation  is  given  by  the  energy  balance 

At/*  =  At/‘*  +  AP  (4.4) 

where  AF*  is  the  critical  energy  release.  It  may  be  assumed  that  there  exists  a  material 
surface  energy  per  unit  area  of  crack  opening  7  ,  in  which  case 

Ar*  =  7A^*  (4.5) 

where  7  is  the  analogue  of  the  critical  energy  release  rate  which  is  usually  defined  as 
the  derivative  of  the  energy  release  with  respect  to  crack  area.  It  can  be  shown  from 
(4.1)-(4.4)  that  the  critical  energy  release  can  be  expressed  entirely  in  terms  of  internal 
stresses  (see  HasHIN  [22]): 

Ar  =  l  f  AaSAtrdV  (4.6) 

V 

For  following  applications  this  criterion  can  be  transformed  to  integration  over  cracks 
surface.  We  can  rewrite  it 

AF*  =  -  J  A<rSA<T  dV  =^J  Aa  AedV 

^  V  V 

=  ^  J  AT  An  dS  (4.7) 

s 

Note,  that  according  to  the  assumption  that  the  traction  on  existing  cracks  remains  zero, 
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the  integration  must  be  evaluated  only  over  surfaces  of  newly  developed  cracks.  Thus 

AP  =  -  [  AT  An  dS  (4.8) 

^  «/ 

ASc 

or 

AP  =  -  /  (t  -  T°)  (u  -  u°)  dS  =  -^  J  T°  (u  -  u°)  dS  (4.9) 

The  surface  integral  for  each  crack  is  confined  to  the  two  congruent  crack  surfaces  and 
thus  each  component  of  the  normal  appears  twice,  with  opposite  signs  (n'*'  =  — n  ,  see 

Figure  4.1). 

Therefore 

Ar  =  -i  /  T"  (u  -  u“)  dS=-\  J  T“+  (u+-u“)  dS 

ASc  ASt 

J 

5c" 

=  -  J  T°(u--u+)  =  ^  J  T°[u]dl/  (4.10) 

ASt 

where  [u]  is  the  displacement  jump  across  the  crack  surface, known  as  crack  opening 
displacement  (COD),  and  the  integration  must  be  evaluated  once.  By  means  of  index 


Figure  4.1:  Crack  faces. 
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notation  (4.10)  may  be  rewritten  in  the  following  form 


7A>U  =  i  /  M  rfS  (4.11) 

ASc 

By  this  is  meant  that  the  energy  release  can  be  expressed  in  terms  of  stresses  in  an 
undamaged  material  in  the  craeks’  locations  and  displacements  of  the  crack’  faces  (COD), 
of  course,  after  development  of  the  craeks.  It  should  be  noted  here  that  kinetic  energy 
cLSsociated  with  the  dynamics  of  crack  growth  has  been  neglected. 

4.2  Fiber  Reinforced  Materials  under  Transverse  Load 

Generally  speaking,  for  a  statistically  homogeneous  material  body,  location  of  new  cracks 
and  magnitude  of  applied  forces  or  temperature  which  lead  to  crack  appearance  depend 
on  loading  conditions.  For  example,  small  concentrated  forces  or  local  heating  may 
produce  microcracks  in  bounded  areas  around  them.  In  order  to  define  the  critical 
loading  conditions  as  composite  material  properties,  we  assume  here  that  the  material 
remains  statistically  homogeneous  after  formation  of  the  cracks  . 

Consider  an  infinite  unidirectional  fiber  reinforced  composite  material  body  subjected 
to  a  homogeneous  temperature  change  </?  and  an  external  applied  homogeneous  traction 
Ti{S)  =  aijTij  in  the  plane  normal  to  the  fibers.  For  effective  transverse  isotropy,  o-.j 
may  be  chosen  in  the  principal  coordinate  system  dn  =  <ti,  0^22  =  o'a?  ^12  =  0  which, 
if  required,  can  be  separated  to  hydrostatic  part  p  =  (di  +  d2)  /2  and  deviatoric  part 
T  =  (di  —  d2)  /2.  Then  the  tensile  radial  stresses  in  the  interface  between  the  matrix 
and  the  fibers  may  produce  interfacial  debonding  around  certain  fibers.  The  further 
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Figure  4.2:  Fiber  reinforced  material  under  transverse  load. 

consideration  will  be  restricted  to  the  case  of  perfect  debond,  when  there  is  no  contact 
area  between  phases  after  development  of  the  interfacial  cracks,  and  the  fiber  surface 
is  free  from  traction.  For  cylindrical  surface  of  fibers  it  is  convenient  to  consider  the 
problem  in  polar  coordinates.  Then  only  the  radial  component  of  n  does  not  vanish  on 
the  interface  (see  Fig.4.2),  and  (4.11)  can  be  rewritten  in  the  form: 

7A/l.  =  i  /  (4.12) 

“aSc 

where  m  denotes  the  matrix  and  /  denotes  the  fiber.  Thus  we  need  to  solve  two  problems: 
first,  we  have  to  find  the  radial  cr®,.  and  shear  stresses  in  the  interface  between  the 
matrix  and  fibers  before  these  fibers  will  separate  from  the  matrix,  and  second,  the 
displacements  of  the  matrix  and  the  separated  fibers  surfaces. 
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4.3  An  Isolated  Fiber  in  an  Infinite  Matrix 


For  one  isolated  fiber  or  dilute  concentration  of  fibers  when  the  interaction  between  fibers 
is  neghgible,  the  interfacial  stress  and  strain  fields  may  be  described  precisely.  Consider 
first  the  case  of  perfect  bond.  The  general  two-dimensional  plane  strain  solution  in  polar 
coordinates  for  transverse  isotropic  elastic  phases  assumes  the  following  forms. 

1.  For  isotropic  load  p  and  temperature  input  (p; 


xir  =  (B  +  C/p'^)  r; 

=  2kB  -  2GCIp^  -  2ka<p 


(4.13) 


where  k  is  the  transverse  bulk  modulus,  G  is  the  transverse  shear  modulus,  a  is 
the  thermal  expansion  coeflficient,  r  =  r fa  where  a  is  the  radius  of  the  fiber  and  B, 
C  are  constants  which  may  be  found  from  the  following  boundary  conditions: 


The  boundary  condition  on  the  external  surface  oo)  =  p 

The  displacement  is  continuous  on  the  interface  (a)  =  (a)  (4-14) 

The  traction  is  continuous  on  the  interface  (a)  =  (a) 

The  additional  condition  C/  =  0  is  needed  to  avoid  singularity  at  r  =  0. 
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2.  For  shear  component  r: 


Ur  =  {Ap^  +  B  +  Cjp"^  +  jD/p^)  r  sin  20; 

ue  =  {^-^Ap^  A  B  A  ^  cos  20; 

1  +  ^  (4.15) 

<7,,  =  2GiB-  I P^  -  Ip*)  sin  20; 

OrB  —  2(j  A  B  A  (7/ p^  -{•  ZD  j p*)  cos  20. 

Here  (  =  Gfk  ,  A,  B,  C  and  D  are  constants,  which  again  may  be  found  from  the 
boimdary  conditions  at  infinity  and  continuity  at  r  =  a  : 

—*  oo)  =  r  sin  20; 
cr^™^(r  — ♦  oo)  =  r  cos  20; 
u(’”>(a,0)  =  4^)(a,0); 

(4.16a) 

4’”>(a,0)  =  4’”)(a,0); 

4-)(a,0)=orW(a,0); 

and  Cf  =  Df  =  0  to  avoid  singularity  at  r  =  0. 

Note  here  that  according  to  the  well  known  result  of  Eshelby,  the  stress  field  in  an  isolated 
fiber  embedded  in  an  infinite  body  imder  homogeneous  traction  is  also  homogeneous, 
and  thus  we  can  obtain  from  (4.15d)  that  Aj  —  0,  and  the  six  equations  (4.16a)  may 
be  reduced  to  five.  The  same  may  be  directly  obtained  from  solution  of  (4.16a).  The 
solution  of  three  linear  equations  (4.14)  for  Bf,  Am-,  Bm  and  any  five  linear  equations 
from  (4.16a)  for  jB/,  Am,  Bm,  Cm  and  Dm  yields  the  following  expressions  for  interfacial 
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stresses  and  displacements: 


=  €rP  +  (llrT  sin  20  + 
o’rfl  =  qlgTcos29-, 

“r  =  {Kp  +  K'’'  2^  + 

—  IgT  a  cos  2^; 


(4.17) 


where  q^r-,  qrri  4  a-nd  Ig  are  the  interfacial  stresses  cr^r,  <^rB  and  the  interface  displacements 
Ur,  ug  respectively  due  to  a  unit  loading  input,  the  nature  of  which  is  noted  by  superscript 
index:  the  hydrostatic  pressure  p,  the  shear  stress  r  and  the  temperature  tp.  These 
coefficients  are  functions  of  the  matrix  and  fiber  properties  only: 


<lrT=9 


1  + 

P  +  6’ 


Prr 


^rr 


=  €b  =  2p 


1  +  6 


2(7mP 


1  +p(l  +  2^m)’ 

«/). 


P  +  6 


(4.18) 


and 


fp  _  6  1  +  ;t  _  /T  _  _ 1  + 

^  2(?mP  +  6’  ^  <Jm  1  +p(l  +  2^m)’ 


(4.19) 


where  g  =  GfJGm- 

After  the  crack  development  cTrr  and  <Trg  vanish  on  the  interface  and  this  situation 
may  be  simply  obtained  from  (4.18)  by  substituting  Gf  —  kf  =  0,  thus  assuming  a  void 
instead  of  a  fiber.  Analogously  the  displacements  of  the  matrix  surface  may  be  obtained 
from  (4.19)  by  putting  Gf  =  kf  =  0  .  This  leads  to  the  following  expressions  for  the 
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matrix  surface  displacements  (4.19c, d): 


“r  =  i^rP  +  sin  2^  +  lr^)a', 
Ug  =  /Jt  a  cos  20; 


with  new  coefficients: 


=  +U;  =  =  +U 


l’?  =  OCm 


(4.20) 


The  fiber  surface  freely  expands  due  to  temperature  change,  so  =  0  and  —  ajacp. 
It  results  in  the  crack  opening  displacement^: 


[ttr]  =  {l^P  +  sin  20  +  {a,n  -  ocf)(p}  a; 
[u^]  = /Jra  cos  20. 


(4.21) 


Now  replacing  p  and  t  by  (<ri+o-2)/2and  (<7i  -<T2)/2  respectively  and  substituting  (4.17) 
and  (4.21)  into  the  debond  criterion  (4.12)  with  help  of  AAc  =  2TvaH  and  dSc  =  aHdO, 
where  H  is  the  length  of  the  fiber,  one  can  obtain: 


^(^Icr  +  ^2cr)  +  “  <^/)^7^cr  +  2C'(ricr<7'2cr 

+2D  (a™  -  OCf){aicT  +  0'2cr)7’cr  = 


(4.22) 


^Generally  speaking,  the  separated  fiber  can  move  freely  in  the  formed  pseudovoid.  But  this  has  no 
influence  on  the  inner  elastic  energy  and  hence  may  be  omitted  from  consideration. 
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where  A,  B,  C,  D  depend  on  the  elastic  properties  of  the  matrix  and  fiber: 


A  = 


^(1  + 


2G„ 


+ 


[4(5 +  ^/)  1  +5'(1  + 2^m)J  ’ 


B  =  2G^ 

n( 

C 


5(1  + 


2Gm 
_  9  i'^+U) 
2  (^  +  ^/) 


[4(5^  +  ^/)  1  +  fl^(l  +  2^m) J  ’ 


(4.23) 


The  formula  (4.22)  relates  the  critical  external  forces  and  temperature  change  which 
produce  interface  debond,  to  phases  properties,  the  critical  energy  release  rate  7  and  fiber 
radius.  The  paxameter  7  is  most  probably  controlled  by  chemical  bond  effectiveness,  the 
fiber  surface  treatment  and  process  of  composite  fabrication. 

In  spite  of  the  fact  that  interfacial  stresses  do  not  depend  on  the  fiber  radius,  the 
critical  loads  and  associated  critical  stresses  in  the  interface  decrease  with  increasing 
of  the  fiber  size,  which  means  that  fibers  of  larger  radius  axe  more  sensitive  to  crack 
formation.  Moreover,  it  follows  that  for  a  given  external  mechanical  load  and  temperature 
there  is  a  critical  fiber  size  above  which  cracks  will  first  appear  and  it  can  be  calculated 
from  (4.22).  Note  here  that  these  conclusions  are  in  conflict  with  the  concept  of  ultimate 
stress  which  is  a  property  of  material,  but  not  of  geometry.  According  to  the  present 
approach,  although  the  magnitudes  of  the  stresses  are  independent  of  fiber  size,  the  total 
stress  energy  depends  on  the  volume  of  material  under  the  influence  of  these  stresses, 
which  in  turn,  depends  on  the  fiber  size.  That  is,  the  larger  the  fiber  radius,  the  larger 
the  stressed  volume  both  within  and  around  the  fiber,  and  thus,  the  larger  the  stress 
energy  that  is  associated  with  the  fiber  (quadratic  in  the  radius).  The  same  may  be  said 
about  the  stress  energy  change  due  to  debonding.  On  the  other  hand,  the  energy  release 
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is  linearly  proportional  to  the  crack  area,  and  consequently  to  the  fiber  radius.  It  follows 
from  here  that  for  a  larger  fiber  size,  smaller  stresses  are  required  for  the  stress  energy 
change  to  reach  its  critical  value  and  for  a  crack  to  appear. 

4.4  Arbitrary  Concentration  of  Fibers 

The  above  result  leads  us  to  conclude  that  in  the  case  of  arbitrary  concentration  of 
fibers,  pseudovoids  may  form  in  a  gradual  manner  aiccording  to  the  size  distribution  of 
fibers.  After  the  first  groups  of  fibers  produce  craeks  in  the  interface  around  them  under 
increasing  load  and  temperature,  the  phase  constitution  of  the  composite  changes;  the 
third  phase,  pseudovoids,  appears  in  the  places  of  debonded  fibers.  Further  increase  of 
loads  will  cause  new  interfacial  cracks  to  appear  and  hence  the  volume  fraction  of  the 
fibers  will  decrease  and  the  pseudovoids  volume  fraction  will  increase.  This  process  of 
voids  formation  may  continue  until  all  fibers  separate  from  the  matrix,  or  in  other  words, 
the  volume  fraction  of  the  voids  is  equal  to  the  initial  volume  fraction  of  the  fibers.  Now 
the  composite  again  contains  two  phases:  the  matrix  and  the  pseudovoids  instead  of  the 
fibers. 

We  restrict  the  following  consideration  to  the  particular  case  when  all  fibers  have 
the  same  radius  a  and  they  separate  from  the  matrix  under  the  same  external  load 
and  temperature.  In  order  to  use  the  criterion  (4.12)  the  generalized  self  consistent 
scheme  approximation  (GSCS)  will  be  applied  to  calculate  the  interfacial  stresses  and 
displacements  for  two  extreme  cases:  a  matrix  with  embedded  fibers  of  volume  fraction 
t)2,  and  a  matrix  containing  voids  of  the  same  volume  fraction  ua-  In  spite  of  the  fact  that 
the  interface  displacements  and  stresses  obtained  in  framework  of  GSCS  are  approximate. 
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Figure  4.3:  Generalized  self  consistent  scheme  (GSCS). 

the  energy  change  may  be  described  in  terms  of  changes  of  the  effective  thermoelastic 
moduli,  which  in  turn  agree  very  well  with  experimental  data. 

The  essence  of  the  GSCS  is  that  in  order  to  determine  the  stress  and  strain  fields  in 
the  phases,  the  composite  is  represented  as  a  composite  cylinder,  consisting  of  a  fiber  core 
and  a  matrix  concentric  shell,  embedded  in  an  infinite  homogeneous  material  to  which 
the  effective  properties  are  assigned.  The  radii  of  the  fiber  a  and  the  matrix  shell  b  are 
chosen  in  such  a  manner  that  (a/6)^  =  U2,  which  is  to  say  that  the  volume  fraction  of  the 
fiber  in  the  composite  cylinder  is  equal  to  the  volume  fraction  of  fibers  in  the  composite 
material  (see  Fig.4.3). 

When  this  structure  is  subjected  to  a  homogeneous  temperature  change  (p  and/or 
hydrostatic  tension  p  is  applied  to  the  external  surface,  the  stress  field  in  material  (0) 
and  consequently  on  the  composite  cylinder  surface  is  homogeneous  (<rf9  =  p),  and 
thus,  the  traction  may  be  applied  immediately  to  the  cylinder  surface  (r  =  b).  This  is 
consistent  with  the  fact  that  the  GSCS  approximation  gives  precisely  the  same  result  for 
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the  etfective  bulk  modulus  k*  as  obtained  by  the  composite  cylinder  assemblage  (CCA), 
according  to  which  the  composite  cylinder  described  above  behaves  just  as  a  homogeneous 
transversal  isotropic  cylinder  with  the  plane  bulk  modulus  k*.  The  foregoing  is  incorrect 
in  the  case  of  a  transverse  shear  input  r,  because  then  the  traction  on  the  composite 
cylindric  surface  does  not  remain  homogeneous  under  a  homogeneous  external  traction. 
For  more  details  about  GSCS  and  CCA  see  CHRISTENSEN  [30]. 

Thus,  the  CCA  in  the  case  of  a  hydrostatic  tensile  load  and  temperature,  and  the 
GSCS  for  the  transverse  shear  will  be  used  for  the  required  interfacial  stresses  and  dis¬ 
placements.  Then  the  general  elastic  solutions  (4.13)  and  (4.15)  are  valid  with  new 
boundary  conditions: 

1.  Uncracked  composite  under  hydrostatic  load  and  temperature  (CCA) 


=  p 

«(/)  (a)  =  «(*")  (a)  (4.24) 

o-C)  (a)  =  (g) 

and  for  pseudovoids 

alfib)  =  p 

(4.25) 

(«)  =  0 

Inserting  equations  (4.13)  into  (4.24)  we  obtain  three  linear  equations  for  unknown 
coefficients  Bj,  and  Cm  (again  (7/  =  0  to  avoid  singularity  at  r  =  0).  The 
solution  results  the  following  expression  for  interface  stress 

^0  _  _ ^m) _ 

‘^lkm{Gm  +  kj)  -f  V2kf{Gm  +  ^ 
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^Gmkfjihj  V\  . 

«MG„  +  k,)  +  v^kt(G^  + 


(4.26) 


After  crack  formation,  recalling  that  =  ayayj,  which  means  free  expansion  of 
fiber  due  to  temperature  change,  equations  (4.25)  give  us  the  desired  COD 


[Ur]  = 


{Gm  +  fcm) 
"i-kjnG mV\ 


P+  (ttm 


-aj)ip 


a 


(4.27) 


2.  For  the  undamaged  composite  subjected  to  pure  shear  (GSCS)  the  boundary  con¬ 
ditions  are 

cr^P(r  oo)  =  r  sin  20, 

—*■  oo)  =  r  cos  20, 
uWI(a,«)  =  u<“)(o.«), 

4^*(“.^)  =  4“’(<>.^). 

u<")  (i.  «)  =  «(»>  (6, «),  ^  _ 

(4.28) 

“«“’(t.^)  =  “f*(^^). 

4;'(a,«)  =  <7W(<>,«). 

<^J»*(“.^)  =  4"’(“.^). 

<.(r>(4,«)=<T<?(6,«), 

(*>.«)  =  <'!?(*.«) 

with  the  additional  restriction  Cf  =  Df  =  0  to  avoid  singularity  at  r  =  0.  Compar¬ 
ing  equations  (4.28a, b)  to  (4.15c,d)  it  is  seen  that 


The  insertion  of  (4.15)  into  (4.28)  yields  eight  linear  algebraic  equations  for  the 
eight  unknown  coefficients:  Af,  Bj,  Cm-,  Dm-,  Co  and  Do-  In  addition, 

Christensen  k  Lo  [29]  have  proved  that 

(70  =  0 


This  precise  condition  has  been  found  on  the  basis  of  the  assumption  that  the 
average  strain  energy  density  of  an  embedded  composite  cylinder  is  equal  to  the 
average  strain  energy  density  of  the  composite  material.  Hence,  the  set  of  eight 
linear  equations  (4.28c)-(4.28j)  gives  G*  as  well  as  the  stress  field  in  the  interface 
between  the  fiber  and  the  matrix,  which  can  be  written  in  the  form 


<7®,.  =  sin  2^ 
cr%  =qreT  cos  2^ 


(4.29) 


In  a  similar  manner,  the  COD  may  be  obtained  when  a  crack  formed  in  the  interface. 
In  this  case  the  continuity  conditions  at  r  =  a  (4.28c,d,g,h)  must  be  replaced  with 


4r*(o,«)  =  o 


which  means  that  the  inner  surface  of  the  matrix  shell  is  free  from  traction.  Solution 
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of  these  six  linear  equations  gives 


[ur]  =  a  sin  20 
[u$]  =  Ig  a  cos  26 


(4.30) 


Superposition  of  (4.26),  (4.27),  (4.29)  and  (4.30)  yields  the  following  expressions  for 
the  required  interfacial  stresses  and  COD 


<^rr  =  q^rP  +  sin26  + 

<^re  =  qlgT  cos  2$-, 

(4.31) 

[«r]  =  {l^P  +  l^T  sin  29  +  (a^  -  «/)¥>}  a; 

\ug\  =  IgT  a  cos  20; 

where  q^^,  qf^  and  are  defined  by  (4.26)  and  (4.27),  and  ql^,  q^g,  are  obtained 
from  solution  of  GSCS  equations  as  described  above. 

Substituting  (4.31)  into  criterion  (4.12)  and  evaluating  the  integration  over  all  fibers’ 
surfaces,  one  can  obtain 


+  <^2cr)  +  ^{oCm  -  +  ^Caicra^cr 

+2D  (am  -  OCf){(Ti„  +  Cr2cr)<pcr  =  —  (4.32) 

a 


where  coefficients  A,  B,  C  and  D  are  functions  of  the  elastic  properties  of  phases  and 
the  volume  fraction  of  the  fibers. 


17 


4.5  Numerical  Results  and  Conclusions 

The  criterion  (4.32)  can  be  rewritten  in  the  form 

F{cricT,cr2cT-><PcT)  =  l  (4.33) 

where  F{<T\cry<^2cri  *Pcr)  is  a  quadratic  polynomial  function  of  critical  applied  stresses  and 
temperature,  which  tend  to  perfect  debond  in  the  interfaces 

(4.34) 

(^)  (^)  +  +  ^)  (^) 

and  c  =  ^27/a.  In  (cri,cr2,(p)  space  (4.33)  builds  a  second  order  surface  which  we  will 
call  the  ’’debond  SURFACE”. 

In  the  (cri,cT2)  plane  the  curve  F(<Ticr,o'2cr,0)  =  1  is  an  ellipse  (Fig .4.4b)  with  the 
semimajor  axis  tilted  at  an  angle  7r/4  to  the  axes  cri  and  <72.  The  values  of  the  ellipse 
semiaxes  may  be  evaluated  by  substituting  crier  =  cr2cr  for  the  semimajor  axis  (which  cor¬ 
responds  to  hydrostatic  tension  2pcr  =  o'lcr  =  <r2cr),  and  <Ticr  =  —(^2ct  for  the  semiminor 
axis  (for  pure  shear  2Tcr  =  crier  =  —o'2cr)~ 

Thus,  if  the  point  (cri,cr2)  which  represents  the  current  loading  conditions  is  inside 
the  ellipse,  then  there  are  no  interfacial  cracks  in  the  material.  Otherwise,  if  the  loading 
path  reach  the  ellipse  surfac^e,  the  energy  release  becomes  enough  to  produce  the  interface 
debond. 

For  the  most  important  case,  when  aj  ^  a„i  ■,  kj  and  G/  >  (?m,  the  temper- 
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ature  <p  scarcely  changes  the  values  of  the  semiaxes,  but  it  has  influence  on  the  ellipse 
center  position.  When  the  temperature  increases  (p  >  0,  the  critical  stresses  decrease 
and  the  ellipse  moves  into  the  third  quarter  of  the  {(Ti,cr2)  plane  along  the  line  ctj  = 
Conversely,  the  ellipse  center  moves  into  the  first  quarter  of  the  (cri,<T2)  plane  along  the 
line  ai  =  when  <  0.  In  both  these  cases  the  ellipse  remains  symmetric  around  the 
lines  <7i  =  <J2  and  ai  =  — a2  +  2h ,  where  ui  =  =  h  is  the  center  of  the  ellipse.  The 

three-dimensional  debond  surface  is  shown  schematically  in  Fig.4.4a.  It  is  obviously  seen 
that  the  surface  presents  ^  elliptic  cylinder  tilled  to  the  (o-i,<T2)  plane. 

The  debond  surfaces  (<Ticr,  o'2ct)  for  various  volume  fractions  of  fibers  V2  and  temper¬ 
atures  ip  are  shown  in  Fig.4.5  and  Fig.4.6.  The  constituent  properties  used  are: 

Epoxy  as  matrix  Glass  as  fibers 
Em  =  0.4  10®psi  Ej  =  10.5  10®psi 
i/m  =  0.35  Vf  =  0.20 

=  30  10-®  °F-i  af  =  2.8  10"® 

It  is  interesting  to  note,  that  the  ellipse  moving  along  the  line  <Ti  =  <7-2  resembles  the 
phenomenon  of  kinematic  hardening,  and  the  increasing  of  the  ellipse  semiaxes’  values 
with  fiber  volume  fraction  decreasing  is  similar  to  isotropic  hardening,  which  are  both 
well  known  from  the  theory  of  plasticity.  Of  course,  in  the  present  problem  this  occurs 
for  completely  different  reasons. 

It  should  be  noted  here  that  not  all  points  of  the  ellipses  correspond  to  possible  phys¬ 
ical  situations.  For  example,  hydrostatic  pressure  and  cooling  can  not  produce  interfacial 
cracks,  because  in  this  case  the  radial  stresses  in  the  interface  are  negative.  In  order  to 
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exclude  this  type  of  problems,  we  have  to  require  that 

K]  >  0  (4.35) 

for  all  angles  9,  in  such  a  way  that  there  is  no  contact  area  between  the  matrix  and  the 
fibers  (perfect  debond),  or  in  other  words,  the  cracks’  surfaces  must  be  free  of  tractions. 
Substituting  the  expressions  for  COD  (4.31c, d)  into  this  inequality  we  obtain  the  ranges 
of  critical  stresses  and  temperatures  that  satisfy  (4.35).  The  points  of  the  debond  surfaces 
which  are  inside  these  ranges  axe  shown  in  Fig.4.5  and  Fig.4.6  in  black. 

All  the  above  leads  us  to  the  following  most  important  conclusions; 

1.  The  smaller  the  fiber  radius  that  is  chosen  during  fabrication  of  the  composite,  the 
larger  the  critical  external  applied  loads  and  temperature  which  produce  interface 
debond.  For  a  given  load  and  temperature  input,  the  maximal  allowable  fiber  radius 
can  be  estimated  from  the  criterion  (4.32). 

2.  As  can  be  seen  from  Fig.4.5,  composites  with  larger  volume  fraction  of  fibers  are 
more  susceptible  to  interface  debond. 

3.  Decreasing  the  temperature  does  increase  the  critical  mechanical  loads,  see  Fig.4.6. 
This  means  that  residual  thermal  stresses  can  prevent  interface  separation.  On  the 
other  hand,  in  the  case  of  overcooling,  the  appearance  of  positive  tangential  stresses 
in  the  matrix  may  produce  matrix  cracking,  which  is  desirable  to  be  avoided. 

4.  It  may  be  assumed  that  the  critical  energy  release  rate  7,  which  can  be  found  by 
experimental  approach,  is  a  property  of  the  interface,  but  not  of  the  volume  fraction 


20 


of  fibers  or  fiber  radius.  The  higher  the  strength  of  interface  between  matrix  and 
fibers,  the  larger  the  loads  that  must  be  applied  to  produce  the  interface  failure. 

These  results  may  be  used  in  design  of  fiber  reinforced  composite  materials  to  prevent 
the  interface  separation  and  improve  the  composite’s  overall  strength. 
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CONCLUSION 


The  research  investigations  presented  in  this  report  have  been  concerned  with  the 
effect  of  microcrack  distributions  on  the  thermoelastic  response  of  composites  and  with 
prediction  of  the  development  of  such  microcracks.  The  underlying  philosophy  may  be 
termed  micromechanics  of  damage.  This  implies  that  defects  such  as  cracks  are  explicitly 
recognized  as  internal  surfaces  on  which  the  tractions  vanish.  This  in  contrast  to 
continuum  damage  mechanics  where  damage  is  described  as  some  mathematical  entity 
(not  uniquely  defined  )  which  involves  parameters  to  be  determined  by  suitable 
experiments.  Micromechanics  of  damage  is  mathematically  more  difficult  but  it’s  great 
advantage  is  in  that  it  is  clearly  connected  to  the  physics  and  geometry  of  the  problems 
considered  and  it  can  yield  results  of  universal  validity  as  has  been  demonstrated  in  the 
studies  presented  here.  The  prediction  of  damage  in  the  form  of  crack  densities  is  a 
particularly  important  problem  and  it  is  believed  that  such  efforts  should  be  continued. 


